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Preliminary group classification became prominent as an approach to symmetry analysis of 
■ differential equations due to the paper by Ibragimov, Torrisi and Valenti [J. Math. Phys. 

32(11), 2988-2995] in which partial preliminary group classification of a class of nonlinear 
^ ' wave equations was carried out via the classification of one-dimensional Lie symmetry exten- 

^ I sions related to a fixed finite-dimensional subalgebra of the infinite-dimensional equivalence 

algebra of the class under consideration. In the present paper we implement, up to both 
usual and general point equivalence, the complete group classification of the same class using 
the algebraic method of group classification. This includes the complete preliminary group 
classification of the class and finding Lie symmetry extensions which are not associated with 
subalgebras of the equivalence algebra. The complete preliminary group classification is 
based on listing all inequivalent subalgebras of the whole infinite-dimensional equivalence 
• algebra whose projections are qualified as maximal extensions of the kernel algebra. The set 

"j^ I of admissible point transformations of the class is exhaustively described in terms of the par- 

tition of the class into normalized subclasses. A version of the algebraic method for finding 
the complete equivalence groups of a general class of differential equations is proposed. 

1 Introduction 

o - 

00 . The method of preliminary group classification was first introduced in Ref. [Ij and became well- 

'n^ I known due to Ref.r |18j . In the latter paper, partial preliminary group classification was carried 

' out for the class of equations of the form 

o : 

Utt = f{x,Ux)Uxx + gix,Ux), (1) 



(N 



(N 



where f ^ 0. The essence of the approach applied in pLSj is given by the classification of one- 
dimensional extensions of the kernel algebra with respect to a fixed finite-dimensional subalgebra 
of the infinite-dimensional equivalence algebra of the class studied. 

The symmetry analysis of the same class was continued in a number of papers. The interest 
to such studies is stimulated because equations of the form ([TJ are used as mathematical models 
of different continuous media. They arise, e.g., in the theory of elasticity, in particular, in the 
course of modeling of hyperelastic homogeneous rods |21l . 

Thus, in jl3l the partial preliminary group classification of the class ([I]) with respect to 
one-dimensional subalgebras of an infinite-dimensional subalgebra of the equivalence algebra 
was considered. Second order differential invariants of the equivalence algebra were computed 
in [19]. Another direction of investigation for the class ([TJ was initiated in |17| . Instead of 
equations of the form ([T]), related systems of two equations, where the first derivatives of u play 
the role of the dependent variables, were considered and mapped to the form vt = a{x, v)wx and 
wt = b{x,v)vx. For the class of systems of this form, certain properties were investigated within 
the framework of symmetry analysis, including the computation of the equivalence and kernel 
algebras and the compatibility analysis of the determining equations for Lie symmetries. Upper 
bounds for the dimension of Lie symmetry extensions were established for the two cases which 
arose. This study was completed in [22] via exhaustive group classification of such systems by 
the algebraic method. 
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A comprehensive review of the hterature on group analysis of different classes of (1 + 1)- 
dimensional wave equations was presented in [27]. Some of these classes are contained in the 
class ([1]) or nontrivially intersect it. In particular, the simple subclasses of ([1]) singled out by 
the constraints fx = g = and fx=gx = were considered in [35] and [12], respectively. The 
class ([T]) has also a subclass common with the class of nonlinear wave equations of the general 
form utt = Uxx + F{t,x,u,Ux), whose Lie symmetries were exhaustively investigated by the 
algebraic method in [27] . The intersection obviously consists of equations of the form ([T]) with 
/ = 1. Any equation of the form ([T]) is a potential equation for the wave equation of another 
form, vtt = {f{x,v)vx + g{x,v))x, also called the nonlinear telegraph equation [71 [T5]. 

Following the paper [TH], several classes of differential equations were investigated within the 
framework of preliminary group classification. Given a class of differential equations, this ap- 
proach in its essence rests on computing optimal lists of inequivalent subalgebras of the associated 
equivalence algebra and studying the Lie symmetry extensions induced by these subalgebras. 
While in the majority of papers on this subject, including Ref. [TS], only symmetry extensions by 
means of inequivalent subalgebras of a fixed finite-dimensional subalgebra of a possibly infinite- 
dimensional equivalence algebra are considered, we have shown in |llj that this restriction is in 
fact not necessary. Stated in another way, there is no obstacle in studying extensions induced 
by subalgebras of the whole (infinite-dimensional) equivalence algebra. This is, what we have 
called the complete preliminary group classification in opposite to the various partial prelimi- 
nary group classifications, which were carried out e.g. in [181148]. As an example, in [11] we have 
solved the complete preliminary group classification problem for the class of nonlinear diffusion 
equations of the general form ut = f{x,u)u'^ + g{x, u)uxx- 

Moreover, in case when the class is normalized (at least in the weak sense [401 144] ) the same 
approach gives at once the complete group classification, cf. Section [3l This fact was implicitly 
used in various instances. The most classical examples for this finding are Lie's classifications 
of second order ordinary differential equations [30] and of second order two-dimensional linear 
partial differential equations |29j. For numerous modern examples see, e.g. [3l (26 1 (27 1 (28 1 l43 l (52 1 
and references therein. The technique of group classification explicitly based on the notion 
of normalized classes of differential equations was developed in [39l HOl HI] and then applied to 
different classes of Schrodinger equations, generalized vorticity equations, generalized Korteweg- 
de Vries equations, etc. All the above techniques can be interpreted as particular versions of the 
algebraic method. 

The purpose of the present paper is to systematically carry out the preliminary group clas- 
sification of the class of differential equations ([1]) in a similar fashion as in [Tlj and thereby 
to exhaustively solve the complete group classification problem for this class of nonlinear wave 
equations using the partition into normalized subclasses. The version of the algebraic method 
applied in the present paper differs from the Lahno-Zhdanov approach [3l [261 113 123 ISl] as it 
does not involve the classification of low-dimensional Lie algebras but is rather based on classi- 
fications of all appropriate subalgebras of the corresponding equivalence algebra. Note that the 
consideration is local throughout the paper. 

In order to guarantee nonlinearity of equations of the form ([T|), we explicitly include the 
nonvanishing condition 

ifu^,gu^uj + (0,0) 

into the definition of the class to be studied. The reason why we are only concerned with the 
nonlinear case here is that nonlinear and linear equations of the form ([T|) are not mixed by point 
transformations (cf. Remark [7]) and have quite different Lie symmetry properties. Moreover, 
linear wave equations of the form ([T|) were already well investigated within the framework of 
classical symmetry analysis in [H [37] . 

The further organization of this paper is the following: Theoretical background of point 
transformations in classes of differential equations is reviewed in Section [2j This includes the 
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definitions and properties of a class of differential equations, its subclasses, the set of admissible 
transformations, the usual equivalence group and algebra, different notions of normalized classes 
of differential equations, etc. Section [3] contains a concise description of the group classification 
problem together with a discussion on the theory of preliminary group classification and complete 
group classification with the algebraic method. Group analysis of the class ^ is started in 
Section [J] by studying the structure of the equivalence algebra of ([1]) . The computation of the 
determining equations for admissible transformations and the equivalence group of the class ([T]) 
by the direct method is given in Section [5] and Section [21 respectively. The algebraic method 
for the calculation of equivalence groups is first presented and then applied to the class ([T|) 
in Section [71 Analysis of the determining equations for Lie symmetries of equations from the 
class (HI) is presented in Section [8l It gives the kernel algebra of this class and allows us to prove 
that the major subclass of the class ^ is weakly normalized with respect to the equivalence 
algebra of (HI). The group classification of the complement of the subclass is also carried out. 
Completing the study of admissible transformations by the direct method, in Section [9] we 
partition the class ([1]) into two subclasses, which are respectively normalized and semi-normalized 
with respect to the equivalence group of the entire class In this way we prove that the 
class dl]) is semi-normalized. Both the equivalence algebra and the equivalence group are used 
in Section [10] to classify subalgebras of the equivalence algebra that may be used for preliminary 
group classification. The adjoint action of the equivalence group on the associated algebra is 
computed using push-forwards of vector fields, as it was recently proposed in [11]. The final 
calculations related to the complete group classification in the class dl]) and the corresponding 
list of inequivalent Lie symmetry extensions can be found in Section [Til In Section [121 we briefly 
sum up the results of the present paper and make comparative analysis of partial preliminary 
group classification, complete preliminary group classification and complete group classification 
within the framework of the general algebraic method. 

2 Point transformations in classes of differential equations 

To make this paper self-contained, in this and in the next sections we restate some important 
notions from the theory of group classification. More information on this subject can be found, 
e.g. in [SlllMlllQllll]. 

The central notion underlying the theory of group classification is an appropriate definition 
of a class of (systems of) differential equations. In practice, the structure and properties of a 
class of differential equations determines which methods of group classification (e.g. complete 
vs. preliminary, direct vs. algebraic) are the most adapted for it. In short, the definition of a 
class of differential equations comprises two ingredients. 

The first ingredient is a system of differential equations Cg: L{x,ui^p-^,6(^g-^{x,U(^p-^)) = 0, 
parameterized by the tuple of arbitrary elements 0{x, 'U(p)) = {9^{x, U(^p^), . . . , 9''{x, ^'.(p))), where 
X = {x^,...,x^) is the tuple of independent variables and is the set of all dependent 
variables u = (n^, . . . , n"^) together with all derivatives of u with respect to x up to the order p. 
The symbol 0(g) stands for the set of partial derivatives of 9 of order not greater than q with 
respect to the variables x and n(p) . 

The second ingredient concerns possible values of the tuple of arbitrary elements 9. This tuple 
is required to run through the solution set 5 of a joint system (also denoted by <S) of auxiliary 
differential equations 5(x, 0(q/)(x, tt(p))) = and inequalities ^(x, 0(g/)(x, «(p))) ^ 0, in 
which both x and n(p) play the role of independent variables and S and S are tuples of smooth 
functions depending on x, U(^p) and 0(q')- The nonvanishing conditions T, ^ might be essential to 
guarantee that each element of the class has some common properties with all other elements of 
the same class, such as the same order p or the same linearity or nonlinearity properties. Thereby, 
these inequalities can be the crucial factor in order to solve the given group classification problem 
up to a certain stage. In spite of this, they are often omitted without any reason. 
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Definition 1. The set {Cg \ 9 £ S} denoted by C\s is called a class of differential equations 
defined by the parameterized form of systems Cg and the set S of the arbitrary elements 9. 

An additional problem in defining a class of differential equations is that the correspondence 
9 ^ Co between arbitrary elements and systems (treated not as formal algebraic expressions 
but as real systems of differential equations or manifolds in the jet space J*-^-* which is the space 
of the variables (rE,U(p))) may not be injective. The values 9 and 9 of arbitrary elements are 
called gauge-equivalent if Cq and Cq are the same system of differential equations, i.e., their 
sets of solutions coincide. We formally consider Cg and £g as different representations of the 
same system from C\s- For the correspondence 6 ^ Cg to be one-to-one in the presence of a 
nontrivial gauge equivalence, the set 5 of arbitrary elements should be factorized with respect 
to the gauge equivalence relation via changing the representation for the class £\s- If this is 
not convenient, the gauge equivalence should be carefully taken into account when carrying out 
symmetry analysis of the class £1^ ^20j . 

In the course of group classification of a complicated class of differential equations, it is often 
helpful to consider subclasses of this class. A subclass is singled out from the class C\s by 
attaching additional equations or nonvanishing conditions to the auxiliary system <S. 

Thus, for the class of equations of the general form ([T|) we have the single dependent variable u 
of two independent variables t and x. The associated tuple of arbitrary elements consists of two 
functions / and g whose domains are contained in the related second-order jet space, i.e., in 
the space of t, x and u together with all derivatives of u up to the second order. The indicated 
dependence of / and g only on x and Ux means that the arbitrary elements of this class are 
solutions of the auxiliary system of differential equations 

ft fu fut futt futx fuxx 0' 

9t = gu = gut = gutt = dut^ = du^x = O- 

As we talk about wave equations, we should also impose the inequality f ^ 0. In the present 
paper we study the subclass of equations of the form ([1]) that consists of only really nonlinear 
equations and, therefore, is singled out from the entire class of equations of the general form ([T]) 
by the additional nonvanishing condition {fu^-, QuxUx) 7^ (0,0). It is the set of equations which 
is called the class ([1]) throughout the paper. 

Several properties hold for subclasses of a class C\s- The intersection of a finite number of 
subclasses oi C\s is also a subclass in C\s-, which is defined by the union of the additional auxiliary 
systems associated with the intersecting sets. At the same time, the complement C\si = C\g, 
of the subclass £1^' in the class £1^ is also a subclass of C\s only in special cases, e.g., if the 
additional system of equations or the additional system of nonvanishing conditions is empty (cf. 
Remark [2]) . Namely, if the subset S' of arbitrary elements is singled out from S by the system 
5[ = 0, . . . , 5^, = then the additional auxiliary condition for S' is |5j^p -|- • • • -|- l-S^/p 7^ 0. If 
S' is defined by the inequalities S'^^ 7^ 0, . . . , S^, 7^ then the additional auxiliary condition for 

5' is • • • s;, = 0. 

A point transformation in a space is an invertible smooth mapping of an open domain in this 
space into the same domain. Given a class C\s of differential equations, point transformations 
related to C\s form different structures. 

Let Cg and Cq be elements of the class £1^. By T{9,9) we denote the set of point transfor- 
mations in the space of the variables (x, u) mapping the system Cg to the system Cg. In this 
notation, the maximal point symmetry (pseudo)group Gg of the system Cg coincides with T{9^ 9). 
If T{6,6) 7^ 0, i.e. the systems Cg and £g are similar with respect to point transformations, 
then T(0, 9) = ip^ o Gg = Gg o if^^ where is a fixed transformation from T{6, 6). 

Definition 2. The set of admissible transformations of the class C\s is given by T = T{C\s) = 
{{9,9,^)\ 9,9eS,^ET{9,9)}. 
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If the number m of dependent variables is equal to one, instead of point transformations one 
can consider more general contact transformations [36] in the same way. 

The notion of admissible transformations \39\ HI] is a formalization of the notion of form- 
preserving \24:\ [25] or allowed [52] transformations. First descriptions of the sets of admissible 
transformations for nontrivial classes of differential equations were given by Kingston and Sopho- 
cleous [2 3 J for a class of generalized Burgers equations and by Winternitz and Gazeau [52j for 
a class of variable-coefficient Korteweg-de Vries equations. An infinitesimal analogue of the no- 
tion of admissible transformations was proposed and studied by Borovskikh in [8]. In terms of 
equivalence groups and normalization properties of subclasses (see Definitions H] and [S]) , the sets 
of admissible transformations were exhaustively described for a number of different classes of 
differential equations which are important for application, such as nonlinear Schrodinger equa- 
tions [43l m] , variable-coefficient diffusion-reaction equations [Ml [50] , generalized Korteweg-de 
Vries equations including variable-coefficient Korteweg-de Vries and modified Korteweg-de Vries 
equations [46j, systems of (1 -|- l)-dimensional second-order evolution equations [45], generalized 
vorticity equations ^40j, etc. 

Definition 3. The (usual) equivalence group = G^{C\s) of the class C\s is the (pseudo)group 
of point transformations in the space of (x, 'U(p), 9) which are projectable to the space of (x, 
for any < p' < p, are consistent with the contact structure on the space of {x, u^p-^) and preserve 
the set S of arbitrary elements. 

Recall that a point transformation (p: z = (p{z) in the space of the variables z = {zi, . . . , z^) is 
called projectable on the space of the variables z' = {zi^ , ■ ■ ■ , Zi^,,), where 1 < h < ■ ■ ■ < ik' < k, 
if the expressions for z' depend only on z'. We denote the restriction of (p to the 2:'-space as 
ip\z'' z' = (p\z'{z'). A point transformation <I> in the space of {x,uijy^,6)^ which is projectable to 
the space of (x,U(p/)) for any Q < p' < p, \s consistent with the contact structure on the space 
of (x,U(p)) if ^\[x,u^p)) is the p-th. order prolongation of 

Each transformation $ from the equivalence group G~ (i.e., an equivalence transformation 
of the class C\s) induces the family of admissible transformations of the form {0 ,^0 ,^\(^x y^) 
parameterized by the arbitrary elements running through the entire set S. Roughly speaking, 
C is the set of admissible transformations which can be applied to any 9 £ S. 

There exist several generalizations of the notion of equivalence group in the literature on 
symmetry analysis of differential equations, in which restrictions for equivalence transformations 
(projectability or locality with respect to arbitrary elements) are weakened [201 l32 l [ii j [ ^ [50]. 

The common part = G'~^{C\s) = CleeS^^ °^ ^ ^ 5, is called the kernel of the 

maximal point symmetry groups of systems from the class C\s [36j. The following folklore 
assertion is true (see, e.g., [TTl [44|). 

Proposition 1. The kernel group G^ of the class £\s is naturally embedded into the (usual) 
equivalence group G~ of this class via trivial (identical) prolongation of the kernel transforma- 
tions to the arbitrary elements. The associated subgroup of G~ is normal. 

Properties of described in Proposition [1] were first noted in different works by Ovsiannikov 
(see, e.g., [38] and [361 Section II. 6. 5]). Another formulation of this proposition was given in [311 
p. 52], Proposition 3.3.9. 

As the study of point transformations of differential equations usually involves cumbersome 
and sophisticated calculations, instead of finite point transformations one may consider their 
infinitesimal counterparts. This leads to a certain linearization of the related problem which 
essentially simplifies the whole consideration. In the framework of the infinitesimal approach, 
a (pseudo)group G of point transformations is replaced by the Lie algebra g of vector fields on 
the same space, which are generators of one-parametric local subgroups of G. 

In particular, the vector fields in the space of (x,u) generating one-parametric subgroups of 
the maximal point symmetry (pseudo) group Ge of the system Ce form a Lie algebra Qg called 
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the maximal Lie invariance algebra of the system Cg. Analogously to symmetry groups, the 
common part = q'^{C\s) = Heesfle 00' ^ G 5, is called the kernel of the maximal Lie 

invariance algebras of systems from the class C\s. It is the Lie algebra associated with the kernel 
group G^. 

The equivalence algebra 0~ is the Lie algebra formed by generators of one-parametric groups 
of equivalence transformations for the class C\s- These generators are vector fields in the space of 
(x, U(p) , 9) which are projectable to the space of (x, U(^pi-j) for any < p' < p and whose projections 
to the space of (x,M(p)) are the p-th order prolongations of the corresponding projections to the 
space of (x, u). 

An infinitesimal analogue of Proposition [T] is the following assertion. 

Corollary 1. The trivial prolongation of the kernel algebra q'~^ to the arbitrary elements is 
an ideal in the equivalence algebra g~. 

By definition, any element of the algebra g'~' formally has the same form as the associated 
element from g^, but in fact is a vector field in the different space augmented with the arbitrary 
elements. 

It is convenient to characterize and estimate transformational properties of classes of differ- 
ential equations in terms of normalization. 

Definition 4. A class of differential equations C\s is normalized if its set of admissible transfor- 
mations is induced by transformations of its equivalence group G~ , meaning that for any triple 
(9, 9, ip) from T{C\s) there exists a transformation $ from G~ such that 9 = ^9 and f = 

Definition 5. A class of differential equations C\s is called semi-normalized if its set of ad- 
missible transformations is induced by transformations from its equivalence group and the 
maximal point symmetry groups of its equations, meaning that for any triple {9, 9, ip) from 
T[C\s) there exist a transformation $ from G~ and a transformation ip from the maximal point 
symmetry group Gq of the system Cq, such that 9 = ^9 and ip = ^\(^x,u) ° ^■ 

In other words, a class of differential equations is semi-normalized if arbitrary similar systems 
from the class are related via transformations from the equivalence group of this class. 

Normalized and semi-normalized classes of differential equations have a number of interesting 
properties which essentially simplify the study of such classes. In particular, if the class L\s is 
normalized in the usual sense, its kernel algebra g^ is an ideal of the maximal Lie invariance 
algebra ge for each 9 & S. In general, this claim is not true even if the class is only semi- 
normalized. See Example 1 in 

The above notion of normalization (resp. semi-normalization) relies on the finite admissible 
transformations. A weaker version of normalization is defined in infinitesimal terms ^40j . 

Definition 6. A class of differential equation C\s is weakly normalized if the union and, there- 
fore, the span of maximal Lie invariance algebras Qg of all systems Cg from the class is contained 
in the projection of the equivalence algebra g" of the class to vector fields in the space of 
independent and dependent variables, i.e. 

U 0e C Pg~ (or {Qe\9eS)c Pg~). 

Here by P we denote the projection operator that acts on vector fields of the general form 
Q = C{x,u)dxi + rj"'{x^u)du'^ + ip'^{x,u,9)dgs in the space of variables x, u and 9 yielding the 
vector fields of the form FQ = C^Xi which are defined on the space of variables x and u. 

It is obvious that any normalized class of differential equations is both semi-normalized and 
weakly normalized. 
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In general, the normalization of a class of differential equations can be checked by computing 
the set of admissible transformations of the class and its equivalence group (e.g. using the direct 
method) and testing whether the condition from Definition U] is satisfied. It is often convenient 
to begin with a normalized superclass and construct a hierarchy of normalized subclasses of the 
superclass or a simple chain of such nested subclasses, which contain the class under consid- 
eration [44 [ H5| I46j . The weak normalization property in turn can be verified by finding the 
equivalence algebra of the class and an inspection of the determining equations for Lie symme- 
tries of systems from the class (see the next section). As the computations related to checking 
weak normalization involve solving of only linear partial differential equations (in contrast to the 
computations using the direct method of finding equivalence and admissible transformations), 
they can be realized in an algorithmic way even for quite cumbersome classes of multidimen- 
sional partial differential equations. At the same time, the established presence of the usual 
normalization property is more useful and allows one to obtain deeper results than involving its 
weak infinitesimal analogue. 

3 Algebraic method of group classification 

Now that we have introduced necessary notions related to point transformations within classes 
of differential equations, we can go on with the general discussion of the framework of group 
classification in some more detail. 

The solution of the group classification problem by Lie-Ovsiannikov for a class C\s of differ- 
ential equations should include the construction of the following elements: 

• the equivalence group G~ of the class £1^, 

• the kernel algebra = q'^{C\s) = Clgi^sQe of the class C\s, i.e., the intersection of the 
maximal Lie invariance algebras of systems from this class, 

• an exhaustive list of G~-equivalent extensions of the kernel algebra g'~' in the class C\s, 
i.e., an exhaustive list of G~-equivalent values of 9 with the corresponding maximal Lie 
invariance algebras qo for which 7^ 0^. 

More precisely, the classification list consists of pairs {S^, {qq, 9 G 5-y}), 7 G F. For each 7 E F 
C\s^ is a subclass of £[5, Qg 7^ g^ for any 9 £ and the structures of the algebras Qg are 
similar for all 9 £ S^. In particular, the algebras Qg, 9 € S^, have the same dimension or display 
the same arbitrariness of algebra parameters in the infinite-dimensional case. Moreover, for any 
9 £ S with Qg g^ there exists 7 G F such that 9 £ mod . All elements from IJ7Gr '^7 
are C-inequivalent. Note that in all examples of group classification presented in the literature 
the set F was finite. 

The procedure of group classification can be supplemented by deriving auxiliary systems 
of differential equations for the arbitrary elements, providing extensions of Lie symmetry, cf. 
Remark [2j In other words, for each 7 G F one should explicitly describe the subset of S 
which is the union of G~-orbits of elements from S^. Although this step is usually neglected, it 
may lead to nontrivial results (see, e.g., [9]). 

If the class C\s is not semi-normalized, the classification list may include equations similar 
with respect to point transformations which do not belong to The knowledge of such 

additional equivalences allows one to substantially simplify the further symmetry analysis of 
the class C\s- Their construction can be considered as one further step of the algorithm of 
group classification [20l [42l |50] . Often it can be implemented using empiric tools, e.g., the fact 
that similar equations have similar maximal invariance algebras. A more systematical way is to 
describe the complete set of admissible transformations. 

In practice, the procedure of group classification within the Lie-Ovsiannikov approach can 
be realized by implementing a few consecutive steps. 
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Given a class £[5, it is convenient to start the procedure by the computation of the equivalence 
algebra. This can be done either using the infinitesimal method [H [36] or simply by deriving 
the set of generators for the one-parametric groups of the equivalence group, provided that 
the latter is known. Computing the equivalence algebra independently from the equivalence 
group is important, as it gives a test and a tool for the calculation of the equivalence group. In 
particular, often only the connected component of unity in the equivalence group is found using 
the knowledge of the equivalence algebra. The equivalence algebra also plays a distinct role in 
the course of applying the algebraic method of group classification. 

The most powerful tool for the construction of the equivalence group, which is the next 
step of the procedure, is the direct method involving finite point transformations. Such a 
construction can be understood as the final stage in the preliminary investigation of the set 
of admissible transformations of the class C\s and allows finding both continuous and discrete 
equivalence transformations. Due to involving finite point transforms the related calculations 
are cumbersome and lead to a nonlinear system of partial differential equations. An alternative 
approach in order to at least restrict the form of point equivalence transformations is based 
on the condition that any point equivalence transformation induces an automorphism of the 
equivalence algebra, cf. Section [71 

The system of determining equations on the coefficients of Lie symmetry operators of a 
system Cg from the class C\s follows from the infinitesimal invariance criterion [6l |3ll 136] . 
stating that 

Q(p)L(x,n(p),6'(g)(x,n(p)))|^p = 

holds for any operator Q = ^^{x,u)dxi + r]'^{x,u)du^ from gg, where the arbitrary elements 6 
play the role of parameters. In what follows we assume the summation for repeated indices. 
The indices i and a run from 1 to n and from 1 to m, respectively. denotes the standard 
p-th prolongation of the operator Q, 

0<|o|^p 

Di = di + u'^idu^ is the operator of total differentiation with respect to the variable Xj. The 
tuple Q = («!, . . . is a multi-index, ai S NU{0}, |a|: = qi + - • • + The variable of the 
jet space J^p^ is identified with the derivative . . . (9x"", and li^ ^ := du^/dxi. Some 

determining equations do not involve the arbitrary elements and thus can be integrated imme- 
diately. The remaining determining equations explicitly depending on the arbitrary elements 
are referred to as the classifying equations. 

Varying the arbitrary elements 9, we can split the determining equations with respect to 
different derivatives of 9. The additional splitting results in equations for those symmetries 
that are admitted for any value of the arbitrary elements and form the kernel of maximal Lie 
invariance algebras. 

The further analysis of the determining equations is usually much more intricate. The classi- 
fying equations are inspected for specific values of the arbitrary elements 9, which give extensions 
of the solution sets of the determining equations, associated with symmetry extensions of the 
kernel algebra. The sets of values found for 9 should be factorized with respect to the equiva- 
lence relation requested. Still, it is the complexity of this analysis that led to the development 
of a great variety of specialized techniques of group classification, which are conventionally 
partitioned into two approaches. 

The first method is the direct compatibility analysis and integration of the determining 
equations up to the equivalence relation, that depends on the values of the arbitrary elements. 
It is mostly suitable for classes with arbitrary elements that are constants or functions of single 
arguments. Algorithms of group classification that are realized in present day's computer algebra 
packages for the calculation of Lie symmetries are based on this method [21 [101 El ED [53] • 
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The other method is of algebraic nature. It is based on the following two properties: For each 
fixed value of the arbitrary elements the solution space of the determining equations is associated 
with a Lie algebra of vector fields. Additionally, if systems of differential equations are similar 
with respect to a point transformation then its push-forward relates the corresponding maximal 
Lie invariance algebras. This is why any version of the algebraic method of group classification 
existing in the literature involves, in some way, the classification of algebras of vector fields up 
to certain equivalence induced by point transformations. The key question is what set of vector 
fields should be classified and what kind of equivalence should be used. 

It is obvious that for each equation Cq from the class C\s its maximal Lie invariance algebra 
09 is contained in the union = UeecSflf- definition of implies that this set consists 
of vector fields for which the system of determining equations is consistent with respect to 
the arbitrary elements with the auxiliary system of the class C\s- Therefore, the set can 
be obtained at the onset of group classification, independently from deriving the maximal Lie 
invariance algebras of equations from the class £1^. As it is not convenient to select linear 
subspaces in the set in the general case, we can extend to its linear span g^^ = (g^l^ € S), 
but fortunately we often have g^ = g^^ Via push-forwarding of vector fields, equivalence (resp. 
admissible) point transformations for the class C\s induce an equivalence relation on algebras 
contained in g^. Such an algebra is called appropriate if it is the maximal Lie invariance algebra 
of an equation from the class C\s- We should classify, up to the above equivalence relation, 
only appropriate algebras. They satisfy additional constraints. The simplest restriction for 
appropriate subalgebras is that each of them contains the kernel algebra g^. The condition that 
the algebras are really maximal Lie invariance algebras for equations from the class is more 
nontrivial to verify. 

Substituting the basis elements of each appropriate algebra obtained in the course of the 
algebra classification into the determining equations gives a compatible system for values of the 
arbitrary elements associated with Lie symmetry extensions within the class C\s- Solving the 
last system completes the group classification within the most general framework of the algebraic 
method. This whole construction is based on the following assertion: 

Proposition 2. Let Si he the subset of S that consists of all arbitrary elements for which 
the corresponding equations from C\s are invariant with respect to the same algebra of vector 
fields, i = 1,2. Then the algebras q^{C\s-^) and q^{C\s2) o-re similar with respect to push- 
forwards of vector fields by transformations from G~ (resp. point transformations) if and only 
if the subsets Si and S2 are mapped to each other by transformations from G~ (resp. point 
transformations) . 

If the class £\s is weakly normalized, the union g^ (resp. the span g^^) is well agreed with 
equivalence. As a result, the algebraic method is appropriate for complete group classification of 
the class C\s. This is not the case when the main part of g^ does not lie in the projection Pg'". 
Then the approach of preliminary group classification [H [18] is relevant to give a partial solution 
of the group classification problem for the class >C|5 by the algebraic method. Preliminary group 
classification essentially rests on the following two propositions (they were first formulated in [18] 
in the particular case of the class ([T]); see |llj for their general formulation and proofs): 

Proposition 3. Let a be a subalgebra of the equivalence algebra g~ of the class C\s, o. C g"", and 
let 6^{x,U(^r)) & S be a value of the tuple of arbitrary elements 6 for which the algebraic equation 
9 = 0'^(x,ii(r)) is invariant with respect to a. Then the differential equation Cgo is invariant with 
respect to the projection of a to the space of variables {x,u). 

Proposition 4. Let Si be the subset of S that consists of tuples of arbitrary elements for which 
the corresponding algebraic equations are invariant with respect to the same subalgebra of the 
equivalence algebra g"" and let be the maximal subalgebra of g~ for which Si satisfies this 
property, i = 1,2. Then the subalgebras ai and 02 are equivalent with respect to the adjoint action 
of G~ if and only if the subsets Si and S2 are mapped to each other by transformations from C . 
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Roughly speaking, in the course of preHminary group classification of the class C\s we classify 
subalgebras of 0~ instead of algebras of vector fields contained in g'-'. Then the objects to be 
classified (subalgebras of g~) are well agreed with the equivalence relation used (C-equivalence). 
If a proper subalgebra s of is fixed and then only subalgebras of s instead of the entire 
algebra g'" are classified up to the internal equivalence relation of subalgebras in s and used 
within the framework of the algebraic method, we call this approach partial preliminary group 
classification. 

In view of Definition [6] and Proposition [3] the following assertion is obvious. 

Corollary 2. For a class of differential equations that is weakly normalized, complete prelimi- 
nary group classification and complete group classification coincide. 

In fact, only certain subalgebras of the equivalence algebra g~ whose projections are contained 
in g^ n Pg'" should be classified. 

Definition 7. Within the framework of preliminary group classification, we will call a subal- 
gebra of the equivalence algebra g~ appropriate if its projection Pa to the space of equation 
variables is a maximal Lie invariance algebra of an equation from the class C\s. 

Appropriate subalgebras of g"" satisfy restrictions similar to those for appropriate algebras 
contained in g^ . As the kernel is included in the maximal Lie invariance algebra of any equation 
from the class, in view of Corollary [T] any appropriate subalgebra o of g~ should contain, as 
an ideal, the trivial prolongation g'^ of the kernel algebra g'^ to the arbitrary elements. The 
condition that the projection Po of a is a Lie invariance algebra of a system from C\s can 
be checked by two obviously equivalent ways: It is sufficient to prove that there exists a value 
9^{x, it(r)) G 5 of the tuple of arbitrary elements 9 for which the algebraic equation 6 = 0^{x, u^^)) 
is invariant with respect to a. The other way is to study the consistence of the system DE^ with 
the auxiliary system of the class C\s with respect to the arbitrary elements. By DE^ we denote 
the system obtained by the substitution of the coefficients of each basis element of Pa into the 
determining equations of the class C\s- Simultaneously we should verify the condition if the 
projection Pa is really the maximal Lie invariance algebra for some systems from jC\s. 

Remark 1. Often the equivalence algebra can be represented as a semi-direct sum of the ideal 
associated with the kernel algebra and a certain subalgebra. To obtain preliminary group classi- 
fication in this case, we in fact need to classify only inequivalent subalgebras of the complement 
of the kernel ideal. Projections of these subalgebras to the space of equation variables will give 
all possible inequivalent extensions of the kernel. This was the case for the class of generalized 
diffusion equations investigated in pJ]. In the present paper, the situation will be different, see 
Remark [9l 

The importance of semi-normalization of a class of differential equations for the optimal 
solution of the group classification problem for this class is connected with the following property 
of semi- normalized classes. 

Proposition 5. For a class of differential equations that is semi-normalized, the group classi- 
fication up to equivalence generated by the corresponding equivalence group coincides with the 
group classification up to general point equivalence. 

In other words, in a semi-normalized class of differential equations there are no additional 
equivalence transformations between cases of Lie symmetry extensions which are inequivalent 
with respect to the corresponding equivalence group. This results in a clear representation of the 
final classification list. As normalized classes of differential equations are both semi-normalized 
and weakly normalized, it is especially convenient to carry out group classification in such classes 
by the algebraic method. This is why the normalization property can be used as a criterion for 
selecting classes of differential equations to be classified or for splitting of such classes into 
subclasses which are appropriate for group classification. 
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4 Equivalence algebra 



The equivalence algebra of the entire class of equation of the general form ([T]) was already 
computed in [18] • It coincides with the equivalence algebra of the class considered in the present 
paper, which consists of purely nonlinear equations of the above form. This is why we only 
represent generating elements of this algebra in a convenient form and refer the reader to [18] 
for more details. The equivalence algebra g~ of class ([T|) is generated by the vector fields 

= udu + u,du^+gdg, V = tdt-2fdf-2gdg, V' = dt, 
V{lp) = ipdx - (fxUxdu^ + 2ip^fdf + ifxxUxfdg, (2) 
= i'du + i^xdu, - ipxxfdg, = tdu, T"^ = t^du + 2dg, 

where (p = ip{x) and tp = ip{x) run through the set of smooth functions of x. The nonvanishing 
commutation relations between the these vector fields are exhausted by 

[^(V'),P"] = a(V), [T\V^]=T\ \T\V^\ = T\ 
[p*, j-i] = j-i, [P*, 7-2] = 27-2, 

In fact, in ([2]) we present only the projections of generating elements of 0~ to the space 
of {t,x,u,Ux, f, g) instead of the whole elements, which are vector fields in the space of 
(t, X, ti(2), /, 5)- The terms of generating vector fields, which are associated with derivatives 
of u, can be computed via prolongation from the coefficients of dt, dx and du and, therefore, 
are not essential. However, it is necessary to include the terms with du^ in the representation 
of these vector fields in order to ensure proper commutation relations between them. Moreover, 
the derivative Ux is a significant argument of the parameter-functions / and g and hence the 
minimal space on which equivalence transformations can be correctly restricted is the space of 
the variables {t,x,u,Ux, f,g)- This is why at least the projections to the same space should be 
given for vector fields from g^. 

The form (3.16) of the equivalence algebra given in |18j differs from ([2]). Namely, the operators 
^(1) and Gix) were singled out from the family {Giip)}- Additionally, we combined the operators 
from [18j to separate scalings with respect to u, which gives simpler commutation relations 
between generating vector fields. 

In order to compute the complete equivalence group of class ([1]) by the algebraic method (cf. 
Section [7]) we need to know a set of megaideals of the equivalence algebra 0~ of this class. 

Given a Lie algebra g, a megaideal i is a vector subspace of q that is invariant under any 
transformation from the automorphism group Aut(3) of 3 [5^, IT]. That is, we have Ti = i for a 
megaideal i of 3, whenever T is a transformation from Aut(0). Any megaideal of q is an ideal 
and a characteristic ideal of g. Both the improper subalgebras of g (the zero subspace and g 
itself) are megaideals of g. If ii and 12 are megaideals of g then so are ii + 12, ii H 12 and [ii, 12], 
i.e., sums, intersections and Lie products of megaideals are again megaideals. If i2 is a megaideal 
of ii and ii is a megaideal of g then i2 is a megaideal of g, i.e., megaideals of megaideals are 
also megaideals. The centralizer of a megaideal is a megaideal. All elements of the derived, 
upper and lower central series of a Lie algebra are its megaideals. In particular, the center 
and the derivative of a Lie algebra are its megaideals. The radical r and nil-radical n (i.e., the 
maximal solvable and nilpotent ideals, respectively) of g as well as different Lie products, sums 
and intersections involving g, r and n ([g,r], [r, r], [g,n], [r, n], [n, n], etc.) are megaideals of g. 

Here we only prove that, roughly speaking, the centralizer of a megaideal is a megaideal. 

Proposition 6. If i is a megaideal of g then the centralizer Cg(i) of i in g is also a megaideal 

ofg. 
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Proof. Consider an arbitrary T e Aut(g), v € i and w G Cg(i). Then [Tif,?;] = [1w,11~^v] = 
1[w,1.~^v] = as T~^t> G i and hence ['w,%~^v] = 0. This means that Tw € Cg(i), i.e., Cg(i) is 
a megaideal of g. □ 

Let = 0~ . It is easy to compute the fohowing megaideals of 0~ : 

g' = {T\V{^),g{i;),T\T'), Q" = {V{^),g{i;),T'), 0'" = (P^, 
Cg(0"') = {V\v\g{l),T\T'), CgK0"') = {V\g{l),T\T'), 

c,,{5") = m),J'\J''), v = (g(i),j-i), z,, = {g{i)), 
Rg = {v'',v\v\g{ij),j^\T^), Rg., = {gm, 

where a', Ra, Zq and Ca(b) denote the derivative, the radical and the center of a Lie algebra a 
and the centralizer of a subalgebra b in 0, respectively. We present proofs only for the last two 
equalities. 

The linear span Si = {T>^,T>^,T'^,g{'ip),J-^,J-'^) obviously is a solvable ideal of 0. Moreover, 
it is the maximal solvable ideal of 0. Indeed, suppose that Si C i and i is an ideal of 0. Then 
there exists a smooth function ( x which does not identically vanish such that the vector field 
T>{(^) belongs to i. As i is an ideal of 0, for an arbitrary smooth function ip of x the commutator 
[D^QjV^if)] = 'D{C,Lpx — Cx'-p) belongs to i. If C, is not a constant function, we define the following 
series of operators: 

R^'' = k-^[V{0,V{C''+% R^^ = k-\R^-^^^,KJ-^'^+\ j,k = 1,2,.... 

It is possible to prove by induction that R^-^'^ = V{C^{CC^f''^) / 0, j, A; = 1, 2, . . . . Moreover, 
as R^^ € i, we have R^^ G i*^-'-', j,k = 1,2, ... , i.e., i^-'^ 7^ {0} for any nonnegative integer j. 
This means that the ideal i is not solvable. If ^ is a constant function, we can set C, = 1. We 
choose any smooth function (/? of x with ipxx ^ and denote ipx by C- As the commutator 
['D(l), = T>{C,) belongs to i, the consideration for the previous case again implies that the 
ideal i is not solvable. Therefore, Si is really the maximal solvable ideal of 0, i.e., Rg = Si. 

The linear spans2 = {g{ip)) is an Abelian and, therefore, solvable ideal of g'". The maximality 
of this solvable ideal is proved in the same way as for Si. Hence Rgm = S2. 

The same megaideals can be obtained in different ways. For example, (^(1)) = Zg/ = Zg///. 

To find more megaideals, we can calculate the automorphism group Aut(m) of a finite- 
dimensional megaideal m and then determine megaideals of m as subspaces of m invariant with 
respect to Aut(m). In the course of calculating the automorphisms we can use knowledge about 
simple megaideals of m. Consider m = C0(0'") = (^(1), J'^ P*, P*). Then 

m' = {g{l),F\F\v'), m" = (a(l),J-^), Zn, = (^(l)), C^{xn") = {g{l),F\F^). 

The presence of the above set of nested megaideals is equivalent to that for any automorphism A 
of m, its matrix in the basis {g{l), F^ , F'^ ,T>^} is upper triangular with nonzero diagonal 
elements. In particular, 

AV* = ai4^(l) + 024-7"^ + as^F'^ + 044^^*, 

AV' = 015^(1) + a25F^ + 0357-2 + 0457^* + 055^^*, 

where 044055 ^ 0. As [P*, D*] = "P* and A is an automorphism of m, we should have [AV^, AV^] = 
AV^ . Collecting coefficients of basis elements in the last equality, we derive a system of equations 
with respect to a's which implies, in view of the condition 044 7^ 0, that 055 = 1, 034 = 0, 
024 = 044035 and oi4 = 044025 — 045024. As a result, we find one more megaideal {g{l),F^ ,V^) 
of m and, therefore, 0~. We will use this megaideal in the course of the computation of the 
complete equivalence group of class by the algebraic method in Section [71 
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5 Preliminary study of admissible transformations 



The infinitesimal invariance criterion allows finding of all continuous equivalence transformations 
by means of solving a linear system of partial differential equations. In order to determine the 
complete point equivalence group (including both continuous and discrete equivalence transfor- 
mations) and the set of admissible transformations, it is necessary to apply the direct method. 
We will start our consideration with a preliminary investigation of the set of admissible trans- 
formations, which will give relevant information also on the equivalence group of the class ([T]). 
That is, we directly seek for all point transformations 

i=T{t,x,u), x = X{t,x,u), u = U{t,x,u), (3) 

for which the Jacobian J = d{T, X,U)/d{t,x,u) does not vanish, that map a fixed equation of 
the form ([T]) to an equation of the same form, 

% = /(x, Ux)uxx + g{x, Ux). 

To carry out this transformation in practice, it is necessary to find the transformation rules for 
the various derivatives of u with respect to t and x. In order to obtain them we apply the total 
derivative operators and D^^, respectively, to the expression u(t,x) = U{t,x,u), assuming 
i = T(t, x, u) and x = X{t, x, u). This gives 

%Dir + MsDiX-Di[/ = 0, 
uPxT + UxDxX - T>xU = 0, 

%(DiT)2 + 2ui-{V)tX){V)tT) + Uxxi^tXf + up^T + u^B^X - D^U = 0, 
Uu(PxTf + 2%^(D,X)(D,T) + Uxx(DxXf + ^D^T + u^bIx - = 0. 

Solving the last two equations for uu and Uxx, respectively, and substituting the results into ([I]), 
we obtain 

Uui^tTf + 2u^-{DtT){DtX) + Uxxi^tXf + u^V'T + UxV'X - V'U 
= f{urt{BxTf + 2ui-{DxT){BxX) + Uxxi^xXf + UiV'T + UxV'X - V^U) (4) 
- giufTu + UxXu - Uu), 

where we use the notation V* = du + 2utdtu + u^duu and = dxx + "^Uxdxu + u^duu and 
additionally have to set u^^ = fuxx + 9 wherever it occurs. As the derivative Uf^ does not 
appear in the transformed form of equations from the class ([1]), the associated coefficient in ^ 
vanishes, i.e. 

{Tt + TuUt){Xt + Xuut) = f{Tx + TuUx){Xx + X^Ux). (5) 

Eq. ([5]) involves only original (untilded) variables and is a polynomial in ut- Therefore, we can 
split it with respect to ut by collecting the coefficients of different powers of this derivative. 
(Note that we cannot as directly split Eq. ([5]) with respect to the derivative Ux, which is an 
argument of the function /.) As a result, we derive that 

-2- TuXu = 0, (6) 
TuXt + TtXu = 0, (7) 
TtXt = f{TxXx + (TuXx + TxXu)ux). (8) 

Multiplying Eq. ([7]) by (resp. Xu), we obtain in view of Eq. ([6]) that T^Xt = (resp. 
TfXu = 0). We apply the trick with the multiplication by Tu (resp. Xu) also to Eq. ([8]) and 
take into account the equations TuXu = 0, TuXt = and TtXu = already derived and the 
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inequality / 7^ 0. This gives equations which involve no arbitrary elements and hence can be 
further split with respect to Ux- Therefore, these equations are equivalent to the equations 
T^Xx = and XyT^ = 0, respectively. The system T^Xt = 0, T^X^ = 0, T^X^ = (resp. 
XuTt = 0, XyTx = 0, XuTu = 0) implies that = (resp. Xu = 0) since otherwise the 
Jacobian J of the point transformation ([3]) vanishes. The condition 

Tu = Xu = 

means that any admissible point transformation of the class ([T]) is fiber-preserving. In view of 
this condition, Eqs. ([6]) and ([7]) are identically satisfied and the remainder of Eq. ([8]) is 

TtXt = fTxXx. (9) 

After substituting u^^ = fuxx + g, we can also split (jl]) with respect to Uxx, which gives, in 
view of Tu = Xu = 0, the equation 

//T,2 + = fifX^ + X,^). (10) 

Unfortunately, the direct splitting with respect to other derivatives in Eq. ^ is not possible. 
The remaining part of (jH) therefore is 

gT^ + u^Ttt + UxXtt - (Utt + 2UtuUt + UuuUt) ^^^^ 
+ UxXxx - (Uxx + '2'UxuUx + Uuuul)) + gUu- 

The additional condition to keep in mind is the nondegeneracy of transformations ([3]) , which in 
view of the conditions Tu = Xu = is reduced to the inequality UuiTtXx — TxXt) / and hence 
{TtXx - TxXt) + and \Ju / 0. 



6 Equivalence group 

At this point, we continue the consideration by computing the equivalence group as it is needed 
even for the analysis of the determining equations for coefficients of Lie symmetry operators and 
the exhaustive description of admissible transformations. In the case of equivalence transforma- 
tions, the arbitrary elements / and g should be varied. We can therefore split the equations ([9]), 
(jlOP and (jlip also with respect to the arbitrary elements. Eq. ([9]) and the nondegeneracy con- 
straint TtXx — TxXt 7^ imply that either Xx = Tt = 0, Tx ^ and Xt or Xt = Tx = 0, 
Xxj^O and Tj / 0. 



For Xx = Tt = 0, Eq. 1^ is simphfied to Xf = ffT^. As the expression for the derivative 



in the new variables is Ux = {TxUi — Ux) /Uu, i.e., it does not involve u^, the equality X^ = ffT^ 
can be split into the two equations Tx = and Xt = 0, which contradict the nondegeneracy 
condition. 

Therefore we necessarily have Xt = Tx = and thus T = T(t), X = X{x), where Xx 7^ 
and Tt 7^ 0. Then Eq. (fTOl) is reduced to fT^ = fX^ and the differentiation of this equation 
with respect to t yields 

2W + /, /^-+/*"^- =0. (12) 

Since Eq. HfTI^i holds for all /, we can split it and derive T^ = 0, Uxt = and Utu = 0. Collecting 
coefficients of in Eq. (jlip we moreover find that Uuu = 0. Taking all the constraints derived 
into account, Eq. (fTT]l reads 

gT^ - Utt = f ( ^"^^^ ^^3: - Uxx - 'iUxuUxj + gUu- 
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Differentiating tliis equation witli respect to u and t allows deriving tliat Uxu = and Uut = 0. 

Collecting all the restrictions derived up to now, any equivalence transformation must satisfy 
the following system of differential equations 

To, = TV = Tff = 0, Xj, = Xf = 0, , , 

_ _ _ (13) 
Uuu — Utu — Uxu — Utx — Uut — 0. 

Integrating the above system in view of the nondegeneracy condition J 7^ 0, we proved the 
following theorem: 

Theorem 1. The equivalence group G~ of the class ([T|) consists of the transformations 

2 C2Ux + Ipx 

t = cit + cq, x = ip[x), u = C2U + c^t +cst + ^p{x), Ux = 



(14) 

f = -^J, 5 = ^ C25 H 'PxxJ - Wxxl + 2C4 

ci Cf V '^x 



where cq, . . . , C4 are arbitrary constants satisfying the condition C1C2 7^ and ip and ij) run 
through the set of smooth functions of x, ipx 7^ 0. 

Comparing the equivalence algebra ^ and the equivalence group ()14p the following corollary 
is evident: 

Corollary 3. The class of equations ([1]) admits three independent discrete equivalence trans- 
formations, which are given by {t, x,u, f, g) i-> {—t,x,u,f,g), {t,x,u, f,g) 1-^ {t, —x,u, f, g) and 
(t, X, u, f, g) H> {t, X, -u, f, -g). 

Theorem [T] implies that any transformation T from G~ of the class ([I]) can be represented 
as the composition 

r = D\ci)7\co)n^W{c2)'J\cA)'y\c:,)%ilj), 

where 



y*(co): 


i = 


t + Co, 


X 


= X, 


u 


= u, 


Ux 


— Uxi 


f = 


f, 


9 


= 9, 


23*(ci): 


t = 


Clt, 


X 


= X, 


u 


= u, 


Ux 


— Uxi 


f = 


cff, 


9 


= Ci^9, 


D{^): 


i = 


t, 


X 


= f, 


u 


= u, 


Ux 


= Ux/^Pxi 


f = 




9 


= 9 + '■PxxUxf Ih^, 


D^{C2): 


i = 


t, 


X 


= X, 


u 


= C2U, 


Ux 


= C2Ux, 


f = 


f, 


9 


= C2g, 


S-^Cs): 


i = 


t, 


X 


= X, 


u 


= U + C^t, 


Ux 


— Uxj 


f = 


/, 


9 


= 9, 


3^(C4): 


i = 


t, 


X 


= X, 


u 


= U + Cit^, 


Ux 


— Uxj 


f = 


/, 


9 


= 9 + 2c4, 


S(V'): 


i = 


t, 


X 


= X, 


u 


= U + ll), 


Ux 


= Ux+ 1px, 


f = 


/, 


9 


= 9- Ipxxf, 



are (families of) one-parameter equivalence transformations and the nondegeneracy requires that 
ciC2(px 7^ 0. These transformations are shifts and scalings in t, arbitrary transformations in x, 
scalings of u, gauging transformations of u with square polynomials in t and arbitrary functions 
of X. 



7 Calculation of equivalence group by the algebraic method 

It is well known that any point symmetry transformation T of a differential equation (resp. a 
system of differential equations) C generates an automorphism of the maximal Lie invariance 
algebra of C via push- forwarding of vector fields in the space of equation variables. This condi- 
tion implies constraints for the transformation T which are then taken into account in further 
calculations using the direct method O [16] . The set of transformations found in the way de- 
scribed form the complete point symmetry group of the equation C including both continuous 
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and discrete point transformations. The above algebraic method can be easily extended to the 
framework of equivalence transformations. A basis for this is given by the following simple 
proposition. 

Proposition 7. Let C\s be a class of (systems of) differential equations, G~ and g~ the equiv- 
alence group and the equivalence algebra of this class (of the same type, namely, either usual 
or generalized ones). Any transformation T from G~ generates an automorphism of via 
push-forwarding of vector fields in the space of equation variables, appropriate derivatives and 
arbitrary elements. 

Proof. Consider an arbitrary vector field Q G 0~ . The local one-parameter transformation group 
G = {exp(e(5)} associated with Q is contained in Then the one-parameter transformation 
group G = {Texp(e(5)T~^}, which is similar to G with respect to T and whose generator is 
%:Q, is also contained in This means that the vector field TIQ belongs to 0~. An arbitrary 
push-forward saves the Lie bracket of vector fields, [%Q,%Q'] = T*[Q,Q'] for any Q,Q' € . 
Therefore, 71 is an automorphisms of □ 

Here we compute the usual equivalence group G~ of class ([1]) by the algebraic method. The 
purpose of this computation is dual: to test the results of the previous section and to present 
an example of applying the algebraic method. The group G~ consists of nondegenerate point 
transformations in the joint space of variables t, x and u, the first derivatives ut and Ux and the 
arbitrary elements / and g, which are projectable to the variable space and whose components 
for first derivatives are defined via the first prolongation of their projections to the variable 
space. Thus, the general form of a transformation T from G~ is 

i=T{t,x,u), x = X{t,x,u), u = U{t,x,u), 

Ui = U\t,X,U,Ut,Ux), Ux = V{t,X,U,Ut,Ux), 



where C/* and are determined via T, X and U and the nondegeneracy condition should be 
additionally satisfied. To obtain the constrained form of T, we will act by the push-forward 7^ 
induced by T on the vector fields ([2]) additionally including the terms with du^ and use megaideals 
of the equivalence algebra 0~ of class ([T]) and restrictions on automorphisms of found in 
Section [H Note that the majority of these objects and properties of g~ are related to the 
finite-dimensional megaideal m = Cg(g"') = , ,V^). Megaideals which are sums 

of other megaideals are not essential for the computation since they give weaker constraints than 
their summands. For example, the megaideal q" is the sum of q'" and Zg/ and hence we do not 
use it in the further consideration. It is sufficient to use the following equalities: 



where a's are constants, (111022033044 ^ and ip and if are arbitrary smooth functions of x. 
The constants o's completed with 055 = 1 and o,j = 0, 1 ^ i < j ^ 5, form a matrix of 
an automorphism of the megaideal m. Tildes over vector fields on the right hand sides of the 



/ = F{t,x,u,ut,Ux,f,g), g = G{t,x,u,ut,Ux, f, g) 



r.J-l = 012^(1) + 022-^\ 

= oi3^(l) + a23T^ + 033-^2, 
T,V^ = oi4^(l) + 024-F^ + 044P*, 

= 015^(1) + 025J^^ + 0357-2 + 045^* + V\ 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



16 



above equations mean that these vector fields are written in terms of the transformed variables. 
The parameter- functions ■ip'^, 'ip'^ and (p'^ are smooth functions of x whose values depend on 
the values oi ip or ip which is indicated by the corresponding superscripts. We will derive 
constraints for T, consequently equating coefficients of vector fields in conditions (fT5]) - (pT]) and 
taking into account constraints obtained in previous steps. As the coefficients of vector fields 
and components of the transformation T associated with the derivatives ut and Ux are defined 
via first-order prolongation involving the similar values related to the variables t, x and u, the 
coefficients of dut and du^ give no essentially new equations in comparison with the coefficients 
of dt, dx and du- This is why we will not equate the coefficients of dut and du^. To have the 
same representation of the final result as in Theorem [H we will re-denote certain values in an 
appropriate way. 

Thus, Eq. (fT5]) implies that Tu = = 0, Uu = C2 and = Gu = 0, where the nonvanishing 
constant an is re-denoted by C2. Then we derive from Eq. ([16]) that tUu = 022?" + 012, i-e. 
T = cit + Co where ci = 02/022 7^ and cq = —012/022, and = Gut = 0. The consequence 
t'^Uu = a33T2 + 023T+ai3 of Eq. (HZD gives only relations between a's. In particular, 033 — C2/c^. 
The other consequences of Eq. pT|) then are Fg = and Gg = C2/cf . The essential consequences 
of Eq. (fT8|) are exhausted by Xt = 0, Ut = a24T + 014 and Ft = Gt = 0. Therefore, X = (p{x) 
and U = C2 + c^t^ + c^t + ipix), where ipx 7^ 0, C4 = a24Ci/2 and C3 = 014 + 02400. 

As we have already derived the precise expressions for the components of T corresponding 
to the variables (cf. Eq. (jl4p ). at this point we can interrupt the computation of equivalence 
transformations by the algebraic method and calculate the expressions for F and G by the 
direct method. At the same time, all the determining equations for transformations from the 
equivalence group G^ of class ([T]) follow from restrictions for automorphisms of the equivalence 
algebra g'". This is not a common situation when the algebraic method is applied. Usually it 
gives only a part of the determining equations simplifying the subsequent application of the direct 
method. See, e.g., the computations of the complete point symmetry groups of the barotropic 
vorticity equation and the quasi-geostrophic two-layer model in [5, Section 3] and ^ Section 4], 
respectively. This is why we complete the consideration of the equivalence group G^ within the 
framework of the algebraic method. 

From Eq. (fT9]) we obtain in particular that tUt = a^^T'^ + 025?" + ois, fFj = F and fGf + 
gGg = G — 035, where 035 = 2c4/cf in view of the first of these consequences. 

Eq. (f20l) implies the equations 

ipUu = i^'^, i^xFu.^=0, i^xGu^-i^xccfGg = i>txF. (22) 

The first and second equations of (j22p are equivalent to = C2ip and F^^ = 0. Then we can 

express the derivative ipf- via derivatives of ip, ipf- = C2(Px^ {^x'ipxx — ^xxi'x), substitute the 
expression into the third equation of ([22]) and split with respect to the derivatives ipx and ipxx, 
as the function ip is arbitrary. As a result, we obtain F = c^'^ifx'^f and Gu^ = C2fx^'PxxF, i.e., 
Gu^ = C2Ci'^(Px^(pxxf ■ The expression for F coincides with the transformation component for / 
presented in Theorem [T] 

The last essential equation for G is given by Eq. (|2ip . Collecting coefficients of dx, du and 
dg in Eq. (I2ip . we have that (p^{x) = (px'P, ^'^(x) = ipx0 and 

ipGx - 0xUxGu^ + 2ipxfGf + ifxxUxfGg = (ptx'^xF - tp'ixF, (23) 

respectively. We proceed in a way analogous to the previous step. Namely, we express the 
derivatives (p"^- and ^p^~ via derivatives of if, substitute the expressions into Eq. ([23]) and split 
with respect to derivatives of p because the function tp is arbitrary. Equating the coefficients of 
(fix leads to the equation fGf = c^"^ ip'^ {c2Ux + ipx 

The simultaneous integration of all the equations obtained for G precisely results in the 
transformation component for g from Theorem [1] 
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8 Determining equations for Lie symmetries 

Suppose that the vector field Q = T{t, x, u)dt + ^(t, u)dx + rj{t, x, u)du belongs to the maximal 
Lie invariance algebra g™'^^ of an equation C: L = from the class ([1]), i.e. it is the generator of a 
one-parameter Lie symmetry group of the equation C The criterion for infinitesimal invariance 
of C with respect to Q is implemented using the second prolongation of Q, which reads 

Qi2) =Q + V'du, + V'^du^ + v''du^ + rf-du,^ + rf-du^.^. 

The coefficients rf" ^ r/**, -rf^ in Q{2) can be determined from the general prolongation formula 
for vector fields, see, e.g. [U |3ll ES]- Using the second prolongation of Q, the infinitesimal 
invariance criterion reads Q(2)-^|l=o = 0, where the notation \l=q means that the condition 
Q(2)^ is required to hold only on equations of class ([1]). Applying the infinitesimal invariance 
condition to the class ([1]) then yields 

"tf^ -[ifx^rffu^xx- frf" -igx-Tfg^^={^ for Utt = fu^x + 9, (24) 

where 

■q^ = D^(7/ - TUt -CUx)+ TUtx + ^Uxx, 
7?^^ =Dl{r]- TUt -£,Ux)+ TUtxx + CUxxx, 
7?" = D?(?7 - TUt -^Ux)+ TUttt + iuttx, 

Dt and D^: denote the operators of total differentiation with respect to t and x, respectively, 
which in the present case of one dependent variable are given by 

T>t = dt + utdu + uttdut + utxdu^ H , 

= 9^ + Uxdu + utxdut + ^Xxdux ~l~ ' ' ' • 

Expanding the infinitesimal invariance condition ()24p we obtain 
D^ri - utB^t^ - n^D^e - 2uttBtT - 2utx,Bt^ 

= f{Blr] - utBlT - u,dU - 2utxB,T - 2u,,B,() (25) 

+ + (D^r] - UtD^T - UxDxi)fuJUxx + ^Qx + (Da;?? - UtDxT - UxDxOdu^ 

where we have to substitute utt = fu^x + g- The above equation can be split with respect to the 
derivatives of utx-, Uxx and ut- Collecting the coefficients of utxUt, utx and UxxUt, we produce 

(.u = 0, Ct = f{rt + TuUx), 2fTu = {Tx + TuUx)fu^. 

Supposing that S,t = 0, it immediately follows from the second equation that r„ = 0. Otherwise, 
for 7^ we can solve the second equation for / and substitute the obtained expression into 
the third equation. After simplification we have that ^tr^ = 0, i.e. Tu = 0. Therefore, we always 
have 

iu = 0, r„ = 0, it = fTx, Txfu^ = 0. 

Further splitting of Eq. (I25p and taking into account the above restrictions gives 
Ut ■ 'nuu = 

Uxx ■ 2{Tt - ix)f + Ux + {rix + {riu - ix)Ux)fu^ = 

Ut ■ 2r]tu - Ttt +Txxf + TxQu^ = (26) 

Rest : 7]tt - ittUx - (Vxx + {2r]ux - Cxx)ux)f + (ilu - 2Tt)g 
- igx - {rix + {r]u - &x)ux)9ua, = 0. 
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The equations = 0, r„ = and r/„„ = neither involve the arbitrary elements f or g nor their 
derivatives. This is why they can be immediately integrated and give restricting conditions on Lie 
symmetries valid for all equations of the form ([TJ. In particular, we have rj = r]^{t, x)u + rf{t^ x). 

In order to derive the kernel of Lie symmetries of class ([T]), we can further split the classi- 
fying part of the determining equations (j26p with respect to the arbitrary elements and their 
derivatives. This immediately gives that the kernel algebra is 

q"" = {dudu,tdu), (27) 

which is a realization of the three-dimensional (nilpotent) Heisenberg algebra. Consequently, the 
Lie symmetries admitted by each equation from the class ([1]) are exhausted by transformations 
of the form (t, x,u) H- (t + ei, x, u + £2 + £3^), where ei, 62 and £3 are arbitrary constants. 

Up to this point the nonlinearity of the equations under consideration was of no impor- 
tance. Only the general form (JT]) was essential. Now we should start to exploit the nonlinearity 
condition (/m^,(7u^u^) 7^ (0,0), which is included in the definition of class ([1]). 

First assume that fu.^ = and therefore Qu^ux 7^ 0. Differentiating the third equation of 
system (|26p with respect to Ux, we then immediately find that = 0. In view of the equation 

= we also have = 0. Upon differentiating the second equation of ([26]) with respect to t 
we obtain r^^ = 0. The third equation of ()26p then implies r]tu = 0. Finally, we differentiate the 
last equation in (p6]) with respect to u and Ux (resp. t and Ux, resp. t). This gives rjxu = (resp. 
r]xt = 0, resp. i^ttt = 0). 

Now we assume that fu^ 7^ 0. In this case, the equation = f^x can be split to yield 
it = Tx =0. The third equation of system ([26]) then implies 2r]tu = tu- Differentiating the 
second equation of (j26p with respect to u we obtain rjxu = 0. The differentiation of the last 
equation in p6p with respect to u then yields r]ttu = 0. In view of the equation 2ritu = tu we 
obviously have Tut = 0. Differentiating the second equation of (j26p twice with respect to t leads 
to rjttx = 0. 

Collecting the results from the above two cases, for the class ([T]), whose definition includes 
the condition {fu^, gu^^u^) (0,0), we always have 

Tu = Tx = Cu= (.t = Vuu = rjxu = Vttx = Tut = 0, 2r]tu = Tu- (28) 

Hence only the second and fourth equations of ()26p essentially involve arbitrary elements and 
are really classifying determining equations for the class ([l]). They must be solved up to the 
equivalence relation induced by transformations from C . Also note that for /^^ 7^ and tu = 
we find by differentiating both these equations with respect to t that ijtx = and ijut = 0. 
This completes the proof of the following proposition: 

Proposition 8. For each equation from class any symmetry operator Q with tu = lies 
in the projection of the equivalence algebra q'" to the space of equation variables, i.e. Q G Pfl'". 

It thus remains to investigate the case in which /^^ 7^ and the corresponding maximal 
Lie invariance algebra g™'^^ = Q^^^{f,g) contains a vector field Q with tu 7^ 0. In view of 
system (j28p the general form of vector fields from g™^^ is 

Q = (02*^ + ait + ao)dt + C{x)dx + ((02* + bi)u + ??°(t, x))9„, 

where the constants oq, ai, 02 and 61 and the functions ^ = ^(x) and r]^ = •rf{t,x), where 
Vttx — 0, are additionally constrained in such a way that the coefficients r = 02^^ + ait + oq, 
^ = £,{x) and rj = {t + bi)u + r]^(t,x) also satisfy the second and last equations of system ([26]) . 
For convenience, we will mark the values of coefficients and parameters corresponding to the 
vector field Q by breve. As 02 7^ 0, by scaling of Q we can set 02 = 1. As the vector field dt 
belongs to g^, it necessarily is in g™^^ for any / and g. Therefore, the algebra g™*^^ contains 
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also the commutator Linearly combining Q with dt and Q = [dt^Q] we can also set 

0-0 = 0,1 = 0. Hence Q = 2tdt + {u + fjt )du- 

Substituting the coefficients of Q into the second equation of system ([26|) . we obtain the 
equation 4/ + {ux + = 0, where /3 = fj^^ should be a smooth function depending at most 

on X since r]^^^ = for any operator Q from g™'^^. The general solution of this equation is 
/ = a{x){ux + /3(x))~^, where a is an arbitrary function of x. Using transformations from 
the equivalence group G~ we can simplify / and set a = ±1 and /3 = 0. If we plug the form 
/ = iu"^ into the second equation of system ([26ll . we obtain that Tt — S,x = "^iVu — ix) + '^Vx^x^ 
for an arbitrary operator from g™^^. From this condition, we can immediately conclude that 
r]x = and ^x = 26i, i.e. ^ = 2bix + 6o for some constant bo. 

As = 0, the substitution of the coefficients of Q into the last equation of system ([26|) gives 
the equation Uxgu^ + 3(7 = with separated variables. Both the sides of this equation are 
equal to a constant which can be set to zero by a transformation 3^(04) from G~. The equation 
Uxdua: +3(7 = is equivalent to the representation g = n{x)u~^, where ji is an arbitrary function 
of X. For this expression of g the last equation of system (I26p takes the form 

- 2binu~^ - {2bix + bo)nxUx^ = 

and the subsequent splitting with respect to Ux implies that r]^^ = and {^{2bix + 6o))x = 0. 
We now distinguish the following cases for values of 60 and 61 depending on a value of /i: 

0. /i is arbitrary. In this case 61 = 60 = 0. 

1. is a nonzero constant. Then 60 is arbitrary and 61 = 0. Using an equivalence transfor- 
mation, we can scale ^ to one. 

2. n = vx^^ mod G~, where is a nonzero constant. (A constant summand of x can be set 
equal to by a shift oi x.) For this value of /i we have 60 = and bi is arbitrary. 

3. yu = 0. This implies that bi and 62 are arbitrary. 

We denote by K. the subclass of equations from the class ([T|, which are G~-equivalent to 
equations with / = ±?x~^ and g = ^{x)u~'^ and by K the complement of this subclass in the 
class ([1]). The above consideration shows that only equations from the subclass K, admit Lie 
symmetry operators that are not contained in Pg". In other words, the following theorem is 
true: 

Theorem 2. The subclass IC of class ([TJ that is singled out by the condition 

if, 9) + {-^'^x^ ,l^{x)u~^) mod 
where [i{x) is an arbitrary function of x is weakly normalized. 

Remark 2. The sets JC and IC of equations are really subclasses of the class ([T]) since the con- 
dition {f,g) = {±u~^ , iJ,{x)u~^) mod G~ and its negation are equivalent to systems of equations 
and/or inequalities with respect to the arbitrary elements / and g. Indeed, by acting on the ar- 
bitrary elements / = iu~^ and g = iJ,{x)u~^ with transformations from G~ and eliminating the 
involved group parameters and the parameter-function fi, we arrive at a system of differential 
equations in / and g characterizing the subclass JC. Namely, the subclass IC is singled out from 
the class ([T| by the system 

Vu^ = i, Wxu.iv\^-w^^iv\u. = o, Wu^uAv\.-WuAv\^u. = o, 

where V = — 4///^^ and W = V^{g + fVx + fxV/2). This implies that the subclass IC as the 
complement of IC is defined by the inequality 
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The above Cases 0-3 represent the complete group classification of equations from the sub- 
class fC up to G~-equivalence. Recall that by the definition of the subclass /C any equation from 
this subclass is G~-equivalent to an equation with / = ±u~'^ and g = ii{x)u~^ . 

Lemma 1. A complete list of -inequivalent Lie symmetry extensions for equations of the 
general form 

utt = ±u~^Uxx + n{x)u~'^, (29) 

where /j, runs through the set of smooth functions depending on x, is exhausted by the following 
cases: 

0. arbitrary ji: 0^ = 0'^ + {t^dt + tudm 2tdt + udu), 

1- = 0"^"^ = 0? + (5.), . 

2. fi = ux'\ u^O: 0>^^- = 0n + (2x9^+n9„), ^^""^ 

3. fi = 0: Q'^^^ = Q'^ + {d,,2xd,+udu). 

Remark 3. We can use equations of the general form 

Utt = 0{x)u~'^Uxx (31) 

as canonical representatives of elements from the class /C instead of ()29p . Indeed, each equation 
from the subclass ()29p is mapped to an equation from the subclass (j3ip by the transforma- 
tion T){ip), where ifxx i ^(px = and 9{x) = ±{ipx{x))^'^ . (Here and in what follows all it and =F 
are consistent with those from Lemma [TJ) In other words, we construct a point-transformation 
mapping |50] between the subclasses (i29l) and (f3T|) which is generated by a family of equivalence 
transformations parameterized by the arbitrary element /i. Hence, mapping and rearrangement 
of the classification list (I30p lead to the equivalent list based on the canonical representative 
form (EH): 



0. arbitrary 9: = g'^ + {t'^dt + tudu, 2tdt + udu), 

1. e = ±e^-: g"'^^ = g'^ + {2dx + udu), . . 

2. 6 = ±\x\^P, py^O: g"^^^ = g^ + (2x5^ + (p+ l)ua„), ^ ' 

3. e = ±l: g'^-''- = g^ + {dx,2xdx + udu). 

Cases 0, 1, 2\^=±i, 2\iy^±i and 3 of the list (pO|) are mapped to Cases 0, 2jp=_.i, 1, 2|p=,^/(j^zpi) 
and 3 of the list (I32p . respectively. Each of the classification lists has certain advantages. Thus, 
the form (j3ip is more compact than (j29p . At the same time, basis elements of the algebras 
presented in the list (pO|) do not depend, in contrast to Case 2 of (f32]) . on equation parameters. 
The equation associated with Case 1 of the list (I30p does not explicitly involve the independent 
variable x, as opposed to its image given in Case 2|p=_i of the list (p2|) whose value of the 
arbitrary element 6 equals ztx"^. 

Remark 4. Due to Theorem [21 to complete the group classification of the class ([I]) it is enough 
to investigate symmetry extensions induced by subalgebras of the equivalence algebra g'". The 
corresponding Lie symmetry generators satisfy the following simplified determining equations: 

Tu = Tx = Ttt = = 6 = Vuu = Vxu = Vtx = Vtu = Vttt = 0, 

Cfx + {{r]u - ix)Ux + r]x)fux = 2(Ca; " Ti)/, (33) 

(.9x + {{ilu - ix)ux + ilx)gux = {riu - 2Tt)g + {CxxUx - ■rixx)f + m- 

Remark 5. Lemma [1] obviously implies that the entire class ([I]) is not weakly normalized. 
This can also be proved without the study of the subclass structure, by the direct computation 
of the union g^ of the maximal Lie invariance algebras of equations from the class ([T]). The 
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set q'^ consists of vector fields of the form r(f, x, u)dt + ^(t, x, u)dx + rj{t, x, u)du for which the 
whole system of determining equations and the nonvanishing condition [fu^jduj^ux) (OiO) 
are consistent with respect to the functions / = f{x,Ux) and g = g{x,Ux)- The consistency 
condition is the joint system of ()28p and 

VtxiVu - Cx) = VxVtu + S,ritxx, 

VtuiCxx + Vxx) = VtxxiVu - '2Tt + 0> (34) 

vtttivu -2Tt) = mivtu - '^Ttt). 

It is clear that is not contained in the projection Pg~ of the equivalence algebra 0~, which 
is associated with the solution set of the system given in the first row of ([33]) . 



9 Set of admissible transformations 

After we have established the equivalence group of the class ([I]), we can describe the set of 
admissible transformations of this class in terms of its normalized subclasses. Theorem [2] and 
its proof give us hints on feasible ways for the classification of admissible transformations. 

First we assume that f^^ = and therefore Qu^u^ 7^ 0. We differentiate Eq. (fTTI) with respect 
to ui and Ux and take into account Eq. Q. From the obtained equation 

TxXx 

u 

and the inequality Qu^u^ 7^ we can conclude that TxXx = 0. Then Eq. Q also implies TtXt = 0. 

Suppose that Tt = 0. Consequently, in view of the nondegeneracy condition of point trans- 
formations we have T^; 7^ and 7^ and therefore Xx = 0. The expressions of ut and Ux 
via ui and Ux take the form ut = {XfUx — Ut)/Uu and Ux = {TxUf — Ux)/Uu, i.e., the expression 
of Ut (resp. Ux) does not involve u^ (resp. n^). We differentiate Eq. ([TT]) twice with respect to 
Ux and once with respect to u. In view of the supposition Tt = 0, this gives {Uuu/U^)u = 0. 
Then we differentiate Eq. (jlip twice with respect to u^: 

9u.u. = 2/^. (35) 

The subsequent differentiation of Eq. ([5^ with respect to u gives the equation {Uuu/Uu)u = 0, 
which together with {Uuu/U^)u = implies that Uuu = 0. Then Eq. ([35|) is reduced to Qu^ux = 
and therefore leads to a contradiction. 

This is why we necessarily have TtXx 7^ and consequently Xt = Tx = 0. In view of Eq. (jlOp . 
we also obtain the equation / = fXx/Tt from which we can conclude by differentiation with 
respect to t that Ttt = 0. It is also evident that = 0. 

Owing to the restrictions derived so far, it is now possible to split Eq. (fTTI) with respect to ut- 
The coefficient of gives Uuu = and that of ut leads to Utu = 0. The rest of Eq. ([TT]) is 

gTf - Utt = f{uxXxx - Uxx - 'iUxuUx) + gUu- 

This obviously implies that gug.ug. 7^ since 7^ 0. We will successively differentiate the 

above rest with respect to three combinations of variables, (u, u^), (t, uj) and t, which gives 
Uxu = 0, Utx = and Uttt = 0, respectively. 

Summing up, for the components T, X and U of admissible transformations of any equation 
with fu^ = and gu^u^ 7^ within the class ([1]) we derive the same system of determining 
equations as in the case of equivalence transformations, cf. ()13p . Moreover, the conditions 
fu^ = and guxUx are saved by the admissible transformations. In this way, we have 
established the following theorem: 
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Theorem 3. The subclass of class ([TJ which is singled out by the constraints fu^ = and 
9uxUa: 7^ is saved by admissible point transformations within the class ([TJ. This subclass 
is normalized and its equivalence group coincides with the equivalence group G~ of the entire 
class ([1]). 

It now remains to investigate the case fu^ ^ 0. Eq. ([9]) immediately impUes that TtXt = 

r,x, = 0. 

Supposing Tx 7^ 0, we obtain that Xx = 0, Xt and hence Tt = 0. In view of these 
conditions Eq. (fTO|) is reduced to X^ = ffT^. Differentiating the last equation with respect 
to Ux leads to the equation ffu^^x = 0, which is equivalent to the equation Tx = 0, contradicting 
the initial supposition. 

Therefore, we have Tx = 0, / 0, = and Xx and Eq. ([10]) reads 



fT? = fXl (36) 
As the transformation rules for the first derivatives are simplified to 

, _ Ut + UuUt ^ _ Ux + UuUx 

T, ' Xx ' 

we can conclude from Eq. ([36]) that fu^ = if and only if = 0. 

Differentiating Eq. ([55]) with respect to u gives {Uxu + UuuUx)fui = 0, and therefore U^u = 
and Uxu = 0. Differentiating Eq. (j36p with respect to t results in the equation 

1: (^^^%^) + ^) + 2/| = 0- (37) 

Taking into account the simplifications obtained so far, we represent Eq. ([TT]) in the reduced 
form 

~rr2 , Ut + UuUt f (Ux + UuUx TT \ , TT 

gTt H Ttt - Utt - 2UtuUt = f I — Xxx -Uxxj + gUu- 

This last equation can be split with respect to uj, giving the equations 



T^Ttt — 2Utu, 
J-t 

gTi + ^Tu -Uu = f r-y-^- x.. -Uxx]+ gU^. 



(38) 



We now distinguish the two cases Tu = and Tu 7^ 0. 

In the case of Tu = 0, the first of the above equations implies Utu = 0. The corresponding 
form of Eq. ([37]) then leads to Utx = 0. Differentiating the second equation of ([38]) with respect 
to t yields Uut = 0. Collecting all the results for this case implies that the transformation 
belongs to the equivalence group G^. 

We now investigate the case of Tu 7^ 0. Then, we solve the first equation of ([38]) with respect 
to Utu/Uu and plug the resulting expression into Eq. ([37]) . This yields 

Ux \ Uxt \ „Tu 



or. 



2Tt V XxJ Xx ' ^ ^ 



Uxt 27t Ux 1 „ 

Ux + ——-—] fui + ^f 
\ ^x J^tt -^x , 
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The difference of the second and third terms in the bracket can be encapsulated as a function 
of X (or, equivalently, x), i.e. we can write {ux + a{x))fu- + 4/ = 0. This imphes that 

f- /^(^) f39) 
^ (u5 + a(5))4- ^ > 

We now differentiate the second equation of ([38]) with respect to u, which gives 

UutTtt 1 fUuTt 



ttu 



Tt 2\ Tt 



where the second equahty holds upon differentiating the first equation in (j38p with respect to t. 
This implies that 

UtuTtt uJTA =0, 



Tt \Tt/t 

which is equivalent to {UuTt/Ttt)^ = 0. Integrating this equation gives an expression for 
Uu = >cTtt/Tt, where x is a constant. We substitute the expression for Uu into the first equation 
of pSj) to obtain 2TtttTt — 3Tj^ = 0. The general solution of the last equation is 

^ _ ait + gp 

were Oj, i = 0, . . . , 3, are constants with 0102 — ao^s 7^ which are determined up to a common 
nonvanishing multiplier. As Ttt 7^ 0, we moreover have 03 7^ and can assume 03 = 1 due to the 
indeterminacy of the constant multiplier. Then we successively gauge 02, ao and oi to 0, 1 and 
by a shift of t, a scaling of t and a shift of t, respectively. All the above transformations belong 
to the group C. In other words, T = 1/t mod C. Plugging the expression obtained for T into 
the equation Uu = >cTtt/Tt allows deriving that Uu = q/t, where g is a nonzero constant. 
Combining Eq. (|36p with the expression for / established in Eq. ([39]) yields 

^ _ T2 aiX)X^ aix) 



{UuUx + Ux + /3(X)X,)4 (u,. + /3(x))4 ' 

where (3{x) := {U^ + p{X)Xx)/Uu and a{x) := Tt^Xxa{X)/U^. Furthermore, upon usmg 
transformations from the equivalence group we can set /3 = /3 = 0, which consequently 
implies that Ux = 0. By means of equivalence transformations, we can also set (3,(3 ^ 1} ^'^id 
as the multiplier relating a and a is strictly positive, we have that a = a. As the transformation 
of X only depends on x it also follows from T^Xx/U^ = 1 that Xx = const. Upon scaling this 
constant and translations of x, which belong to we can choose X = x. Therefore, Tf /U^ = 1. 
As T = 1/t and thus Uu = q/t, this means that ^ = 1, i.e., Uu = 1/t and hence U = u/t + U^{t) 
and Ux = Ux/t. Here U^ = U^{t) is a smooth function arising after integration with respect to u 
and depending only on x in view of the condition Ux = 0. 
The remaining part of Eq. (Illh can be represented as 

9_ ^ttO ^ttO 



^-2C/0-tC/0 =9, (40) 



where U = u/t + U^{t). The differentiation of Eq. ([iO]) with respect to t yields ga-Ux + 3g + 
t^{tU^^ + 2U^)t = 0. The first two terms do not depend on t and the last summand depends 
only on t. Thus, we can separate variables and set t*^(tC/j*j + 2U^)t = — 3x = const, where the 
factor of —3 was introduced for the sake of convenience. Integration of this equation yields 
tUj/f + 2Uf = x/t^ + X, where x = const. The general solution of this equation is U^ = 
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>c/{2t^) — xt/2 — ai/t + (jQ, where cti,(T2 = const. We also have ga-Ux + 35 = ^, which upon 
integration leads to ^ = jl{x)/u^ + k. Plugging these results into Eq. (|10]) gives 
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Using equivalence transformations, we can put k = k = 0. This is why we have / = bju, 
f = 5/u'^, g = and g = d/w^, where 5 = ±1. That is, the equivalence transformations for 
this case reduce to symmetry transformations. 

Owing to the above computations, we can formulate the following theorem: 

Theorem 4. The subclass K, of the class ([1]), that consists of equations -equivalent to equa- 
tions of the form (j29p . is semi-normalized with respect to Any admissible transformation 
in this subclass is generated by G~ or is represented as a composition of the transformations 
i9i,92,Ti), {92,92,T2) and (^2,^3,73), where 9i = {f,g), 92 = {±u-\fiu-^), 9^ = {f,g) and 
Ti, are equivalence transformations and T2 = 1/t is a symmetry transformation of Cg^. The 
complement K, of K, in the class ([1]) (as well as the complement of K, in the subclass of ([1]) 
singled out by the condition ^ 0) is normalized with respect to . The usual equivalence 
group of the subclass JC coincides with G~ . 

Corollary 4. The entire class (JI]) is semi-normalized. Hence the group classification of the 
class ([T|) up to G^ -equivalence coincides with the group classification of this class up to general 
point equivalence. 



Remark 6. It can be proved using the above consideration that the class (|29p is normalized. 
The equivalence group of this class consists of the transformations of the general form 



ait + ao , , . , , ~ ±W\lnA\u-\-b3t-\-b2 _ /i 
t = —— , x = bit + bo, u = —^ — , ^ = 7-, 

031 + «2 «3t + 0.2 bi 

were a,, i = 0, ... ,3, are constants with A = aia2 — 00^3 7^ which are determined up to 
a common nonvanishing multiplier and bi, i = 0, ... ,3, are arbitrary constants with bi 7^ 0. 
The group G^ can be represented as the product of its two subgroups. The first subgroup is 
the ideal associated with the kernel group of the class (j29p and formed by the transformations 
from Gi with 61 = 1 and b^ = 0. The second subgroup corresponds to the subgroup of G~ 
whose elements save equations of the form (|29p and consists of the transformations from Gi 
with 03 = 1 and 02 = 0. This is why the list presented in Lemma [T] is an exhaustive list of Lie 
symmetry extensions in the class (j29p up to both G^^-equivalence and general point equivalence. 

Remark 7. It follows from the above consideration that the entire class of equations of the 
general form ([T]) is partitioned into three subclasses associated with the additional constraints 
fu^ / 0) /«x = and ^ 0, and = gu^u^ =0, respectively. Equations from different 

subclasses of this partition are not mapped to each other by point transformations. This is the 
main reason why it is natural to separate nonlinear equations of the form ([1]) from linear ones, 
which are well studied and form the last subclass. 

Remark 8. In order to simplify calculations, we could use Theorem 4.4b of Ref. [24j, describing 
form-preserving transformations between (1 + l)-dimensional second-order partial differential 
equations of the quite general form uu = H{t,x,u,Ux,Uxx), where H^^^ ^ 0. This theorem 
directly implies the simplest constraints T^ = Tx = = = for admissible transformations 
of the class ([1]), in view of which the coefficients of any Lie symmetry operator Q = rdt+idx+r^du 
of each equation from the class ([I]) satisfy the determining equations t„ = = ^„ = = 0. 
A partial repetition of computations in the present paper was necessary in order to finding the 
appropriate partition of the class ([1]) into subclasses. 
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10 Classification of inequivalent appropriate subalgebras 



In order to classify subalgebras of the equivalence algebra we need to describe the adjoint 
action of the equivalence group C, which consists of transformations of the form (jl4p . on the 
generating vector fields ([2]) of g~. This adjoint action can be determined by solving the Cauchy 
problem 

dw 

— — = ad v(w) := [v,w], w(0) = wq 
de 

for each pair (v,wo) of generating vector fields of 0~, which is equivalent to computing the 
convergent Lie series [34] 

oo „ 

w(e) = Ad(e^")wo := V ^(ad v)"(wo). 

n=0 

An alternative way is the direct computation of actions of transformations from G~ on elements 
of 0~ via pushforward of vector fields by these transformations [llj. Stated in another way, the 
second method uses the usual transformation rule of vector fields under point transformations. 
As this method properly works for infinite-dimensional Lie algebra, we will pursue it below. 

Employing elementary equivalence transformations (cf. the end of Section E]), we can compute 
the nonidentical adjoint actions using the respective push- forwards. This yields 

3-2 (C4)P* = V' + 2c4 -F^ , Vi (ci ) = , 

where 6 = 9{x) is the inverse of the function 9. It should be stressed that there are more 
nonidentical adjoint actions of transformations from on generating vector fields of than 
listed above, namely those related with actions on the trivial prolongation of the kernel 
algebra g^ to the arbitrary elements, which is an ideal in g", and those involving T^cq) and 
3^l{cs). These adjoint actions, however, do not yield simplifications in the course of classification 
of extensions of the kernel algebra. 

We will only classify appropriate subalgebras of g"". Any appropriate subalgebra s of g~ 
should contain g^ = {V*^ , J-^ , G {1)) . For the class ([1]) we have two specific representations of 5, 
which are given by s = g" (Q^ . . . , Q'^) = (V^T^) + (^(1), Q^,...,Q^), where "-h" denotes 
the direct sum of vector spaces, g'^ is an ideal of s (since it is an ideal of the entire g~) and 
{G{1),Q^ , ■ ■ ■ , Q^) is a subalgebra of s. Q^^ . . . , are basis elements from the complement of 

in s and their projections to the space of equation variables yield a proper Lie symmetry 
extension of g'^ in the class ([T]). 

Remark 9. The double representation of appropriate subalgebras is related with the rep- 
resentation of the whole equivalence algebra g~ in the form g~ = g'~' + g, where and 
g = {T>'^,T>^,'D{ip),Gi^),J^) are an ideal and a subalgebra of g~ but the sum is not direct 
even in the sense of vector spaces since g'^ H g = (^(1)). Unfortunately, the equivalence alge- 
bra g"^ does not possess a representation as a semi-direct sum of the ideal g^ associated with the 
kernel algebra and a certain subalgebra, which additionally complicates the group classification 
of the class ([1]). 

This is why it is necessary to classify only subalgebras of g"" which are contained in g and 
contain {G{1)). The classification should be carried out up to Gj^-equivalence, where Gq is a 
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subgroup of formed by the transformations (|14p with cq = C3 = 0. In fact, we will present 
the classification results in terms of extensions of q'^ excluding G{1) from the corresponding 
bases. 

The determining equations for Lie symmetries of equations from the class ([1]) impose more 
restrictions on appropriate subalgebras. 

Lemma 2. s n {V^, Q{ip),T'^) = s n (P*, J-"^) = {0} for any appropriate subalgebra s. 

Proof. Suppose that an appropriate subalgebra s of g"" contains an operator Q = bV^ + Q{tp) + 
cJ-'^, where at least one of the constants b and c or the derivative ipx of the function ijj = il){x) 
does not vanish. Then the operator VQ is a Lie symmetry operator for an equation from the 
class ([T|). Substituting the coefficients of operator Q into the determining equations ([33]) implies 
the following conditions for the arbitrary elements / and g: 

{bux + ipx)fu^ = 0, {bux + Tpx)gu^ = 2c - ipxxf + bg. 

For both the cases 5 7^ and 7^ it follows that fu^ = and gu^u^ = 0, which contradicts 
the definition of class ([1]). The case 6 = 0, V'x = and c 7^ leads to a contradiction. Therefore, 
any appropriate subalgebra does not contain an operator of the form considered. 

Analogously, an operator + cJ-"^, where c is an arbitrary constant, gives the condition 
/ = 0, which is also inconsistent with the definition of the class ([T]). □ 

Lemma 3. dim (s n ^(i/'), -F^)) ^ 2 /or any appropriate subalgebra 5. 

Proof. Suppose that s is an appropriate subalgebra of and dim (s n {'D{ip),Q{ip), T'^)) ^ 2. 
This means that the subalgebra s contains (at least) two operators = T){if'^) + G{ip^) + CjJ-"^, 
where the functions (^*, i = 1,2, should be linearly independent in view of Lemma [21 In other 
words, the projections PQ* of simultaneously are Lie symmetry operators of an equation from 
the class ([T]). By W we denote the Wronskian of the functions ip^ and ip'^, W = if^iPx — V^V'x- 
VF 7^ as the functions (p^ and ip'^ are linearly independent. 

Plugging the coefficients of PQ* into the first classifying equation from the system (j33p gives 
two equations with respect to / only, 

i^lux - ^Dfu^ - p'fx + 2plf = 0. (41) 

We multiply the equation corresponding to z = 2 by (/?^ and subtract it from the equation for 
i = 1 multiplied by ip'^. Dividing the resulting equation by W, we obtain the ordinary differential 
equation 

{ux + /3)/«. + 2/ = 0, 

where /3 = I3{x) := {^'^tpx — ip^'4>^)/W and the variable x plays the role of a parameter. It is 
possible to set /? = by means of an equivalence transformation, S(— /3). Indeed, this trans- 
formation preserves the form of the operators Q*, only changing the values of the functional 
parameters V*- In particular, it does not affect the linear independency of the functions 93*. 
The integration of the above equation for /3 = yields that / = au~'^ , where a = a{x) is a 
nonvanishing function of x. In view of the derived form of /, splitting of equations (j4ip with 
respect to Ux leads to p^ax = and (/J^ = 0, i.e. V'x = and ax = 0. As Q{1) G s, we can assume 
up to linear combining of elements of s that = 0. The constant a can be scaled to a = ±1 
by an equivalence transformation. 

In a similar manner, consider the last equation from system (I33p . taking into account the 
restrictions set on parameter- functions and the form of /. For each Q*, this classifying equation 
gives an equation with respect to g, 

f'xUxgu^ - p'gx = -iPxx'^Ux^ - 2ci. (42) 
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Again, we multiply the equation corresponding to i = 2 with ip^ and subtract it from the 
equation for i = 1 multiplied by , divide the resulting equation by W and thereby obtain that 
9u^ = + ^^^^x^i where /i^ = /^^(x) := aWx/W and //^ = yu^(x) := 2{cnp^ - C2(p^)/W. 

Integration with respect to Ux directly gives g = ^J''^u~^ + In \ux\ + fjP, where fi^ = ijP{x) is a 
smooth function of x. The parameter-function /i^ can be set equal to zero by the equivalence 
transformation D(C)i where the function ( = ({x) is a solution of the equation a(xx + M^Cx = 0- 
Substituting the derived form of g into equations (I42p and splitting with respect to Ux, we find 
that /i^ = 0, ipxx = 0, if^fJ-x = V^xA*^ + 2ci. Therefore, fi^ is a constant and the functions ip^ 
and can be set to 1 and x, respectively, upon linear combining of Q^. Then, we have 
Ux = 2ci, x^x = 2c2 + /i^, i.e. ci = 0, C2 = —1/2 and /x" is a constant that can be set to zero by 
the equivalence transformation 3^{—fj,^/2). 

Summing up, we have proved that any equation of class ([I]) admitting (at least) two opera- 
tors PQ* is G~-equivalent to an equation of the form 

Utt = ±Ux^Uxx + /"^ In \Ux\, 

where /i^ = const. However, the determining equations ()33p in this case yield rjx = 0, ry„ = r^, 
Cxx = 0, fi^rju = 0, r]tt = ^^{Tt — ^x)- This obviously implies that the number of such operators Q* 
cannot exceed two. □ 

Corollary 5. There are two -inequivalent cases of Lie symmetry extensions in class ([I]) 
involving two linearly independent operators of the form PQ*, where = T>{ip'^) + ^("0*) + CiJ-^ , 

1. uu = ±u-^uxx + 2\u\ux\: g'""^ = g^ + (PP(l),PP(x) -PJ-2), 

2. utt = ±u-^uxx: g'"^- = gn + (PP(l),PP(x),PP* + PP"). 

Proof. For ^ 0, we have that rj^ = Tt = 0, £^xx = and, after scaling of /x^ to two by an 
equivalence transformation, rju = —2(,x- This directly gives the first case. If fj,^ = 0, we obviously 
recover the second case. □ 



Now that we have computed the essential adjoint actions and classified all appropriate sub- 
algebras in Corollary [5] for which dim (s fl {'D(ip),Q(ilj), T'^)^ = 2, we should go on with the 
computation of inequivalent appropriate extensions of g^, which contain at most one linearly 
independent operator of the form T>{ip) + G{Tp) + cJ-"^, where if = (p{x) is a nonvanishing func- 
tion. In view of Lemma [2] it is obvious that the dimension of such extensions cannot be greater 
than three. Here we select candidates for such extensions using only restrictions on appropriate 
subalgebras presented in Lemmas [2] and [3l As there exist specific restrictions for two- and three- 
dimensional extensions, we will make an additional selection of appropriate extensions from the 
set of candidates directly in the course of the construction of invariant equations. 

The result of the classification is formulated in the subsequent lemmas. 

Lemma 4. A complete list of -inequivalent appropriate one- dimensional extensions of g'"' 
in g~ is given by 

{V'-G{x)), {V{l) + eF^), ^^^^ 
where e G {0, 1} and a is an arbitrary constant. 

Proof. The classification of the appropriate one-dimensional extensions can be carried out effec- 
tively by simplifying a general element of the linear span {V'^,V^,'D{ip),Q{'>p),J^), 

Q = aiP" + asP* + V{^) + g(V^) + a^T^, 
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using push- forwards of transformations from G~. For this aim, it is necessary to distinguish 
multiple cases, subject to which of the constants or functions (p, ip are nonzero. We note in 
the beginning that owing to the push- forward D^{9) we can always set ip = = const. 

For ai 7^ we can scale the vector field Q to achieve ai = 1. Using the push- forwards of 
a suitable transformation S(x)) we can set ■0 = 0. The further possibilities for simplification 
depend crucially on the value of 02. For 02 = 1/2, the sum I?" + 02 is invariant under the 
push- forward 3"^ (04) and therefore it is not possible to set 04 = 0. The actions of the push- 
forwards of the transformations T>{x) and D" allow scaling of 03 and 04 to {0, 1}. If 04 = 1, by 
denoting 03 = e we obtain the first case from the list (j43p . 

For 02 / 1/2 we can use the push-forward 9'*(c4) to additionally set 04 = 0, which gives, 
jointly with the case 02 = 1/2 and 04 = 0, the second extension listed, where 02 is denoted 
by -p. 

If ai = and 02 7^ 0, we scale 02 = 1 and can use the push-forward 3"^ (04) to set 04 = 0. 
For a3 7^ 0, we can scale 03 = —1 by means of the action of T)'^{c2) and additionally put ■0 = 
upon using the push- forward of the transformation S(x)- If 03 = 0, we have "^a; 7^ in view of 
Lemma [2] and hence we can use the action of T>^{9) to set ip = —x. This gives the third and the 
fourth case of the list (j43p . respectively. 

In case of ai = 02 = but 03 7^ 0, we can set 03 = 1 and use the push-forward 9*(V') for 
a certain to arrive at •0 = 0. The action of I'"(c2) on the resulting vector field allows us to 
scale the coefficient 04 so that we have 04 € {0, 1}, which yields the fifth element of the above 
list of one-dimensional inequivalent subalgebras. 

In view of Lemma [21 the case ai = 02 = 03 = is not appropriate. □ 

Lemma 5. Up to -equivalence, any appropriate two-dimensional extension ofg'^ in which 
contains at most one linearly independent operator of the form T>{ip) + ^(0) -|- cF"^ , belongs to 
the following list: 

(P" + P(l), + (p-" + P(l), P* + g(e")), 

{aiV" + a2V' + a^V{x) + eog{x)+eiF'^,V{l) + e2F'^), 

where b, ai, 02, 03, eo; £1 o,nd £2 are constants with 6 7^ 0, (01,02) 7^ (0,0), (02, ^3) 7^ (0,0), 
(ai,a3,eo) (0;0, 0) and (ai — 2a2 — 03)52 = 0. Due to scaling of the first basis element and 
G''^ -equivalence we can also assume that one of a's equals 1, (ai — 03)69 = 0, (oi — 202)^1 = 0, 
eo, £1 € {0, 1} and £2 G {—1, 0, 1} or, if eq = 0, £2 G {0, 1}. 

Proof. The general strategy is to take two arbitrary linearly independent elements and 
from the linear span (P", 2?*, P((^), ^(■0), J^^) such that s = g'"' -|- {Q^,Q'^) is a five-dimensional 
subalgebra of satisfying the restriction on elements of the form 2?((/5) + g{ip) + cJ-"^ and 
Lemma [2] and simplify and as much as possible by linear combining of elements of 5 and 
push- forwards of transformations from The proof is split into two parts. 

First, we consider possible extensions not involving operators of the form 'D{if) + Q{'ip) + cJ-'^. 
In view of this additional restriction and Lemma O up to linear combining we can take the 
elements and in the initial form 

Qi = P« + V{ip^) + g(V^) + ci-F^, Q2 ^ + V{ip^) + g{^P^) + C2T^ 

where ip^ 7^ 0. We set p^ = 1, = and ci = using D*(^), S*(x) with suitable functions 9 
and X and 3"^(ci), respectively, i.e. = T>^ -\- 1^(1). As the subalgebra s is closed with respect 
to the Lie bracket of vector fields, we have Q^] = V{ipl) + Q{tpl - tp"^) - 02^"^ G (^(1)) and 
hence (^^ = 0, C2 = and ■0^ — ■0^ = const. Integrating the equations for ip^ and 0^ gives that 
(pp' = b and ■0^ = die^ -\- d^ for some constants b, di and do. The constant do can be always 
set equal to zero by linear combining with the operator ^(1) belonging to g^. The further 
simplification of depends on the value of b. lib ^ 0, the push-forward of S(— di6~^e^') does 
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not change and leads to di = 0. If 6 = 0, the parameter di is nonzero in view of Lemma [2] 
and, therefore, can be scaled to 1 by D*((i^^). As a result, we obtain the first two elements of 
the list dH]). 

Now we investigate the case dim (s fl {T>{lp),Q{iI)),F'^)) = 1. Then basis operators of the 
extension of can be chosen in the form 

where (01,02) 7^ (0,0) and ip^ ^ 0. We set (/?^ = 1 and ip'^ = Q using D*(0) and S*(x) 
with suitably chosen functions 9 and respectively. As s is a Lie algebra, we have that 
[Q^,Q^] = V{ipl.) + ^(V'i;) + (oi — 2a2)c2-F^ = 03(52 + Q{d) for some constants 03 and d. 
Therefore, (p^. = 03, (ai — 2a2 — 03)02 = and V'i; = cq- Up to combing with and ^(1) we 
obtain that (p^ = a^x and V'^ = coa;. Up to C-equivalence we can assume that (ai — 03)09 = 
and (oi — 2a2)ci = 0. Indeed, if oi — 2a2 7^ 0, we can set ci = using 3'^(ci) with an appropriate 
constant ci. To set cq = in the case oi — 03 7^ 0, we simultaneously act by S*(co2;) with an 
appropriate constant cq and linearly combine the operator with T>{1). Using push-forwards 
of scalings of variables and alternating their signs, we can independently scale the constant 
parameters cq, ci and C2 and change sings of ci and, simultaneously, cq and C2. Additionally we 
can multiply the whole vector field by a nonvanishing constant to scale one of nonvanishing 
a's to one. The conditions (02, 03) 7^ (0, 0) and (oi, 03, eo) 7^ (0, 0, 0) follow from LemmaO This 
yields the third case of the list ()44|) and thereby completes the proof of the theorem. □ 

Lemma 6. Up to -equivalence, any appropriate three-dimensional extension of in 0~, 
which contains at most one linearly independent operator of the form D{ip) + Qi^) + cF^ , has 
one of the following forms: 

{V^+PiV{x), V'+p2V{x), V{1) + eF'), 

{v'' + v{x) + dg{x),v'-g{x),v{i)), 

where pi, p2 and d are constants, e £ {0, 1}, (pi,P2) 7^ (l^O) and e{pi — 1) = e{p2 + 2) = 0. 

Proof. In view of Lemma[2l any three-dimensional appropriate extension of g^, which contains at 
most one linearly independent operator of the form T){ip)-'rQ{Tp)+cF'^ , is spanned by vector fields 
Qi = p«+p(^i)+g(^i)+ci7-2, q2 ^ pt+p(^2)+^(^2)^^27-2 and Q'^ = V{ip^)+g{^j^)+C3F^ , 
where ip^ and tp^ are smooth functions of x, Ci are constants, 'y?^,'^^ 7^ and (99^, V'^) 7^ (0,0). 

Using 3~^,(ci), T>^{6) and 9*(x) with suitably chosen functions 9 and x of x and, if C3 7^ 0, 
I'^(cJ^), we set ci = 0, ip^ = 1, ip^ = and C3 = e € {0, 1}. The commutation relations of 
with and are 

[Q3, Q1] = V{^1) + ^(V'^.) + C3-F2 = piQ' + diGil), 

for some constants and di, i = 1,2. These commutation relations imply the conditions ip^. = Pi, 
ipl, = di and (pi — l)e = (p2 + 2)e = 0. Therefore, up to combining with and ^(1) we 
obtain ip^ = piX and ip'^ = diX. Then the commutation relation 

[Q^, Q^] = {d2 + P2di - Pid2)g{x) + C2F'^ = 

yields C2 = and P2di = {pi — l)d2- If Pi 7^ 1, we can set di = using 9*(— (pi — l)^^dix) 
and then the equality p2di = {pi — 1)^2 is reduced to c?2 = 0. Analogously, in the case P2 7^ 
we can set d2 = using 9*{—P2^d2x) and then the equality P2di = {pi — l)d2 is equivalent 
to di = 0. Summing up, we always have di = d2 = mod G~ if {pi,P2) 7^ (1)0). This gives 
the first extension in (ll5]l . Otherwise, pi = 1, p2 = and hence e = and ^2 7^ 0. Setting 
(i2 = — 1 by D"(— ^2 ^), we obtain the second extension in (jl5]) . This completes the proof of the 
theorem. □ 
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11 Result of group classification 



To describe equations whose Lie invariance algebras contain the projection Ps of a certain 
appropriate subalgebra s of 0~ to the variable space, we can use two equivalent ways, which 
lead to the same system of partial differential equations in the arbitrary elements / and g: For 
each basis element Q of s we should either substitute the coefficients of PQ into the last two 
equations of system ()33p or write the condition of invariance of the functions / and g with respect 
to Q. Then we should solve the joint system of the equations derived. Simultaneously we should 
check whether the projection Ps is really the maximal Lie invariance algebra for obtained values 
of the arbitrary elements / and g. 

All the candidates for one-dimensional appropriate extensions listed in Lemma |31 are really 
appropriate. For each representative of the list we have an uncoupled system of two equations 
in / and g, which is easy to be solved. As a result, we obtain the following list of equations from 
class ([T|) that admit one-dimensional Lie symmetry extensions of related to 

LL P" + ip* + + ^2 

1.2. V -pV* + V{e): 

1.3. P*-P(l): 

1.4. V^-g{x): 

1.5. V{l) + eT^: 

where a; = x — eln \ux\, £ E {0, 1} and p is an arbitrary constant. Here and in what follows, in 
each case we present only basis elements the corresponding subalgebra of 0~ which belong to 
the complement of the basis {T'^,g{l),T^} of g^. 

Calculations related to two-dimensional extensions are more complicated. We first present 
the result of the calculations and then give some explanations. 

2.1. P« - V{p), - P(l), p^O,: utt = ±e-2-|u,.|2p(u,., + uu^), 

2.2. V^ + Vix), V^-g{x): uu = ±e'^''- (x^u^^ + ux), 

2.3. P« + P* + P(x), P(l): uu = fMu,,, 

2.4. V'^ + V' + V{x) + gix), V{1): -Ui* = ±m^^ + e-"^ 

2.5. 2P" + P* + 2V{x) + g{x) + V{1) : uu = ie^^-u^^ + 2u^, 

2.6. P" + V{x) + g{x), V{1) + e2^^ £2 G {-1, 0, 1} : uu = ie^^-u^^ + e"- + 2£2X, 

2.7. (2 - + (1 - + (2 - q)V{x) + g{x), V{1), g / 0, 1 : 

2.8. i2 + 2p- q)V'' + {1+p- q)V^ + {l + 2p- q)V{x), P(l), 

2.9. (3 + 2p)P" + + (1 + 2p)V{x), + : 
Uu = ±K\^Puxx + e3\ux\P+^/^ + 2x, £3 G {0, 1}, 

2.10. 2(l+p)P" + (l+p)P* + (l + 2p)P(x)+7-2, V{1): = ±hxf ^^^xx + 21n |u^.|, 

2.11. 2P« + P* + 27'^ ^(1) + J'^. Uu = ±u-^Uxx + 2 In \ur,\ + 2x. 

Nontrivial constraints for constant parameters which are imposed by the maximality condition 
for the corresponding extensions are discussed in detail after Theorem [5l 

Cases 2.1 and 2.2 correspond to the first and second spans from LemmaO respectively. For the 
associated invariant equations to have a simpler form, these span are replaced by the equivalent 
spans (D" — D* — D(l)), where p = —b~^, and {V^ + 'D{x), — g{x)), respectively. Note 
that we always can set a constant multiplier of the arbitrary element / to ±1, e.g., by scaling 
oft. 

The third span from Lemma [5] in fact represents a multiparametric series of candidates 
for appropriate extensions, which is partitioned in the course of the construction of invariant 



Utt = f{uj)u~^Uxx + g{uj) + 2 In \ux 
Utt = \ux\'^P{f{uj)uxx +g{uj)ux), 
Utt = e'^'-^{f{ux)uxx +g{ux)), 
Utt = e^'^''{f{x)uxx+g{x)), 
Utt = hux)uxx + 9{ux) + 2ex), 
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equations into Cases 2.3-2.11. Not all values of series parameters give appropriate extensions. 
Additional constraints for parameters follow from the consistence conditions of the associated 
system in the arbitrary elements, 

fx = 0, ((ai - 03)^^ + eo)fux = 2(a3 - 02)/, 

gx = 2e2, ((ai - 03)11^. + eo)gux = (oi - 202)5 - 26203^; + 2ei, 

with the inequality / 7^ and the requirement that the dimension of extensions should not 
exceed two. 

The above partition is carried out in the following way. 

If oi = 03 = a2, the common value of a's is nonzero and we can set it to be equal 1 by scaling 
the first basis elements of the span. We also have that ei = mod G~ and 82 = 0. Depending on 
either eg = or eg = 1 (which is replaced by the equivalent value eo = ~1) we obtain Cases 2.3 
and 2.4, respectively. 

If oi = 03 7^ 02, scaling the first basis elements of the span allows us to set 03 — 02 = 1. 
The parameter eo should be nonzero since otherwise / = 0. Therefore, eo = 1 mod The 
conditions 02 = 1, 02 = and 02 7^ 0, 1 lead to Cases 2.5, 2.6 and 2.7, respectively. In the last 
case we denote 1 — 02 by q. This value should also be nonzero since otherwise the extension 
dimension is greater than two. 

Let oi 7^ 03. Then eo = mod G~ and by scaling the first basis elements of the span we can 
also set 01 — 03 = 1. Introducing the notation p = 03 — 02 and q = ai — 2a2, we obtain that 
oi = 2 + 2p — g, 02 = 1+p — q and 03 = l + 2p — q. The further partition depends on values of e2, 
q and ei. For e2 = the dimension of extension is not greater than two only if either (7 7^ and 
then ei = mod C (Case 2.8) or g = and ei 7^ and then ei = 1 mod (Case 2.10). The 
condition e2 = 1 implies that q = p + 1/2. If additionally either g = ei = Oor(77^0 (and then 
ei = mod G~), we have Case 2.9. Case 2.11 corresponds to the additional constraints g = 
and ei 7^ (i.e. ei = 1 mod C). 

Consider the candidates for three-dimensional appropriate extensions listed in Lemma [6l 
The compatibility of the associated systems in the arbitrary elements, supplemented with the 
inequality / 7^ 0, implies P1+P2 = ^ and d = for the first and the second span of Lemma [H re- 
spectively. The general solutions of these systems up to G~ -equivalence are {f,g) = {±\ux\'^^,0) 
and {f,g) = (±e^"'',0). This gives the following cases of Lie symmetry extensions: 

3.1. (l+p)P"+pP(x), (l+p)P* + P(x), P(l), P7^ -2,-1,0: uu = ±\ux\^Puxx, 

3.2. V"" + V{x) + g{x), V -g{x), V{1): = ±e2"-tx^^. 

Special cases of Lie symmetry extensions in class ([1]) are presented before this section. More 
precisely, all inequivalent equations whose maximal Lie invariance algebras are not contained 
in the projection of the equivalence algebra 0~ to the variable space are listed in Lemma [TJ 
Equations from class ([T]) which are invariant with respect to two linearly independent operators 
of the form PQ*, where = 'D{ip^) + Q{7p^)+CiJ-'^ , are described in CorollaryO For convenience, 
we collect the derived cases in a single table and formulate the final result of group classification 
in the class ([T]) as a theorem. Recall that within the class ([I]) C-equivalence coincides with the 
general point equivalence, cf. Corollary HI 

Theorem 5. All -inequivalent (resp. point-inequivalent) cases of Lie symmetry extensions 
of the kernel algebra in the class ([T]) are exhausted by cases presented in Table [TJ 

In each case of Table [1] we present only basis elements of the corresponding Lie invariance 
algebra that belong to the complement of the basis {dt,du,tdu} of g^. The spans of g^ and the 
vector fields given in cases 1-6 and 9 of Table [U are the maximal Lie invariance algebra of the 
corresponding equations for the general values of the associated parameter- functions / and g, 
but for certain values of these parameter-functions additional extensions are possible. 
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Table 1: Lie symmetry extensions of the kernel algebra = ((9t, du,tdu) for the class ([T]) 



N 


/ 


g 


Basis of extension 


One-dimensional extensions 


1 


f{x — e In \ux\)ux ^ 


g{x - eln \ux\) + 2ln \ux\ 


tdt + 2ea, + 2(u + t^)du 


2 


f(x — eln \ux\)\ux\^^ 


g{x — eln luxDluxf'^Ux 


—ptdt + edx + udu 


3 






tdt — dx 


4 






tdt — xdu 


5 


Rux) 


g{ux) + 2ex 


dx + et^du 


Two-dimensional extensions 


6 


5ux 


g{x)ux^ 


t^dt + tudu, 2tdt + udu 


7 




ue'^'luxl^Pux, iy(p+l)^S 


pdx - udu, tdt — dx 


8 




i/a;e^""=, u ^ S 


xdx + udu, tdt — xdu 


9 


fiUx) 





dx, tdt + xdx + udu 


10 


s 




dx, tdt + xdx + {u + x)du 


11 




2ux 


dx, tdt + 2xdx + {2u + x + t^)du 


12 




e"''+2£2X, e2G{-l, 1} 


xdx + {u + x)du, dx + £2t^du 


13 




e""- , q / 


dx, (1 - q)tdt + (2 - q)xdx + ((2 - q)u + x)du 


14 




M", *) 


dx, {l+p~q)tdt + {l+2p~q)xdx + (2+2p-q)udu 


15 




e|u,|f+i''2+2x 


dx + t^du, tdt + {l+2p)xdx + {3+2p)udu 


16 




2 In lij; 


dx, {l+p)tdt + il+2p)xdx + {2{l+p)u + t^)du 


17 




2 In \ux \ + 2x 


dx + t'^du, tdt + 2{u + t'^)du 


Three-dimensional extensions 


18 






t^dt + tudu, 2tdt + udu, dx 


19 


5ux* 


UX~^Ux^, V ^ Q 


t^dt + tudu, 2tdt + udu, 2xdx + udu 


20 


5|u,p^p/-2,o 





dx, tdt + xdx + udu, ptdt ~ udu 


21 







dx, tdt + xdx + udu, tdt — xdu 


Four-dimensional extensions 


22 


5ux ^ 





t^dt + tudu, 2tdt + udu, dx, 2xdx + udu 



Here S ^ ±1 mod G~ and e £ {0, 1} mod G~. In Case 15 c = mod G~ if p = -1/2. 
*) g / 0, (p,q) + (-1, -1), (-2, -3) in Case 14. 



In the course of collecting cases of Lie symmetry extensions into Table [H they are properly 
arranged. In particular, Cases 1 and 2 of Corollary [5] are merged with Cases 2.10 and 3.1 into 
Cases 16 and 20 of Table [H respectively. As the value p = — 1 is singular for the basis of 
Case 2.10, the bases of Case 2.10 and Case 2 of Corollary [5] are changed in order to be agreed. 
Case 2.6 with £2 = is not included in Case 12 of Table [1] since it is united with Case 2.7 into 
Case 13 of this table. 

Within the algebraic approach used for group classification of the class ([T]) the construction of 
Lie invariance algebras precedes the construction of associated invariant equations. This is why 
the simplification of the form of bases of Lie symmetry extensions, in a certain sense, dominates 
in Table [TJ The form of invariant equations can be slightly simplified if simultaneous minor 
complication of bases of the corresponding Lie invariance algebras are permitted. In particular, 
multipliers equal to two can be removed from arbitrary elements by equivalence transformations 
or re-denoting the parameter p. 

Note that the unique inequivalent case of Lie symmetry extension for which the corresponding 
Lie invariance algebra is of maximal dimension possible for equations the class ([T|) and equal to 
seven. Case 22, is not associated with a subalgebra of the equivalence algebra . 
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Now we discuss nontrivial constraints for constant parameters which are imposed by the 
maximahty condition for the corresponding extensions. 

The equation utt = e'^^\ux\^^{5uxx + I'Ux) corresponding to Case 7 for general values of 
parameters is linear if p = 0. If p 7^ —1, it is reduced by the transformation i = \p + Ij^^^^t, 
X = e~^^^^'^'^\ u = u to the equation Ufj = \ux\'^^iSuxx + i'x~^Ux) with i) = 6 — v{p + 1), which 
coincides with the equation of Case 19 (resp. 20, resp. 22) if p 7^ —2 and P 7^ (resp. p = —2 
and v ^ resp. p = —2 and v = Q). 

The equation utt = e^"^"^ {6x'^Uxx + i^xux) corresponding to Case 8 is similar with respect to 
the transformation i = t, x = x, u = u + xln\x\ — x to the equation n^^ = e^"* {Suxx + {1^ — S)x~^) 
which coincides with the equation of Case 21 \i v = 5. 

Consider the subclass of the class ([1]) associated with the additional constraints fx = 9x = 0, 
i.e., the class of equations of the general form 

Utt = f{Ux)Uxx + g{ux), (46) 

where {fu^^iQu^u^^) 7^ (0)0). Lie symmetries of these equations were comprehensively described 
in [12 ] using no equivalence relation. The selection of Cases 5|e=o, 9, 10, 11, 13, 14, 16, 18, 20, 21 
and 22 from TablelH which are related to equations from the subclass (|46p . represents the exhaus- 
tive list of Lie symmetry extensions in this subclass up to general point equivalence, where the 
algebra = (5t, tdu, dx) given in Case 5|£=o is the corresponding kernel algebra. Theorems[3] 
and [4] imply that the equivalence group of the subclass (j46l) is a subgroup of the equivalence 
group namely, it consists of transformations of the form (jl4p with ipxx = "^xx =0. As 
all necessary shifts, scalings and sign changes of the derivatives uu^ Uxx and Ux are induced by 
transformations from the equivalence group, the majority of constants parameterizing elements 
of the classification list from [12] can be set to appropriate values (0, 1, ±1 or others). In 
other words, these constants are inessential from the point of view of symmetry analysis. As a 
result, the classification list from [12j is reduced to the above selection of cases from Table [H 
excluding a single case given in [12] as Case X. After simplifications by shifts, scalings and sign 
changes of derivatives induced by transformations from , the value of arbitrary elements for 
this case takes the form / = and g = 6ux^, where 6 = ±1. A value of the coefficient 

of Ux^ in the expression for g is not essential. We set it equal to 6 for convenience. The re- 
lated equation uu = 6u~'^Uxx + Su~^ has the six-dimensional maximal Lie invariance algebra 
{dt, du,tdu, dx, e^dx,tdt + udu) and is reduced by the transformation T: i = t, x = , u = u to 
the equation tt^^ = 6u'^'^Uxx, which corresponds to Case 20|p=_i from Tabled! This is why we put 
(p, q) 7^ (—1, —1) as a parameter constraint for Case 14 from Tabled) The presence of one more 
inequivalent case in the course of the classification of the subclass (l46]l up to G^-equivalence in- 
stead of general point equivalence is explained by the fact that the transformation T does not be- 
long to G2 ■ To complete the discussion of singular values of parameters for Case 14 from Tabled! 
we only note that the values p = —2 and q = —3 directly give Case 18 and for the value q = the 
corresponding equation is reduced by the transformation u = u — 1'^/2 to the equation associated 
with Case 20. The same transformation reduces Case 13 with q set to zero formally to Case 21. 
Any equation from the class ([T|) is a potential equation for the equation of the form 

vtt = {f{x, v)vx + g{x, v))x (47) 

with the same value of the arbitrary elements / and g, where the argument Ux is replaced by v. 
Indeed, Eq. (|17|) possesses two inequivalent characteristics of conservation laws, = 1 and 
= t. The potential systems constructed with the simplest conserved vectors associated with 
these characteristics is 



'wl = vt, wj = f{x,v)vx + g{x,v), (48) 
wl = tvt-v, wf = tf{x,v)vx +tg{x,v). (49) 
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We denote tw^ — iiP' by u. In terms of the dependent variables v, and n, the joint potential 
system ([18]) . (|19]) takes the form Ux = v,ut = w^, = f{x,v)vx + g{x,v) which is a potential 
system for system ()48p . i.e., it is formally a second-level potential system of Eq. (|47|) . Hence 
ti is a second-level potential for this equation. Excluding v and from the last system, we 
obtain Eq. ([T]). In order to derive Eq. (j47l) from Eq. ([T]), we should totally differentiate Eq. ([T|) 
with respect to x and replace Ux by v. As the coefficients of any Lie symmetry operator Q = 
rdt + (,dx + r]du of Eq. ([T]) satisfy the determining equations = S,u = Vuu = Vxu = 0, the 
coefficient of in the prolongation of this operator to v according to the equality v = Ux ''\s 
equal to rjx + {rju — Cx)ux and hence does not depend on u. Therefore, Lie symmetries of Eq. ([1]) 
induce no purely potential symmetries of Eq. (f^7|) . 

We checked cases from Table [U using the package DESOLV |10t [5T] for symbolic calculations 
of Lie symmetries, whenever it was possible. 

12 Conclusion 

Results of this paper and those existing in the literature on symmetry analysis of differen- 
tial equations allow us to comparatively analyze different approaches to group classification 
of differential equations (partial preliminary group classification, complete preliminary group 
classification and complete group classification) within the framework of the algebraic method. 
Given a class C\s of (systems of) differential equations with the equivalence group G~ and the 
equivalence algebra 0~, the application of each of the above approaches involves, in some way, 
classification of certain subalgebras of 0~. The essential point is what subalgebras of g~ should 
be classified and what equivalence relation should be used in the course of the classification. 

In the course of partial preliminary group classification, a proper subalgebra 5 of 0~ is fixed 
and then only subalgebras of s are classified. This approach may be relevant only if the sub- 
algebra 5 is noticeable from the physical or another point of view. Hence the choice of such 
subalgebra should be strongly justified which, unfortunately, is often ignored in the existing 
literature on the subject. Differences in the consideration of the subalgebra s instead of the 
whole algebra g~ are especially significant in the case when 0~ is an infinite-dimensional alge- 
bra whereas s is a finite-dimensional subalgebra. A seeming advantage of replacing g"" by s is 
that in general finite-dimensional algebras are much simpler objects than infinite-dimensional 
ones. At the same time, partial preliminary group classification has a few essential weaknesses 
most of which are related to the following fact: As the fixed subalgebra s of g~ is usually not 
invariant under the adjoint action of the equivalence group this group does not generate a 
well-defined equivalence relation on subalgebras of s. This is a reason why subalgebras of s are 
classified up to the weaker internal equivalence on s, which is induced by the adjoint action of 
the continuous transformation group associated with s, instead of G~-equivalence. 

The exhaustive classification of subalgebras up to the internal equivalence is a cumbersome al- 
gebraic problem, possessing no algorithmic solution even for finite-dimensional algebras. In order 
to simplify it, only one-dimensional subalgebras are usually classified which crucially increases 
incompleteness of results obtained in the framework of partial preliminary group classification. 
Although the number of classification cases remains quite large, many of them are inessential 
from the G'^-equivalence point of view, not to mention the general point equivalence. The 
presence of equivalent cases unnecessarily complicates both the solution of the group classifica- 
tion problem and further applications of classification results, e.g., for the construction of exact 
solutions of systems from the class C\s. 

Complete preliminary group classification of the class £1^ is based on the classification of 
subalgebras of the entire equivalence algebra 0~ up to G~-equivalence. As both the objects, 0~ 
and are directly related to the class C\s and well consistent to each other, this approach looks 
as quite natural. For weakly normalized classes of differential equations, it gives an exhaustive 
classification list. Moreover, complete preliminary group classification always is a necessary step 



35 



for complete group classification within the framework of the algebraic method. It is obvious that 
complete preliminary group classification gives a list which is closer to exhausting all possible Lie 
symmetry extensions than any list obtained via partial preliminary group classification. At the 
same time, due to the usage of G~-equivalence which is stronger than the internal equivalence 
on a subalgebra of the former list can contain even a less number of cases than the latter 
one. For example, 33 cases of one-dimensional extensions of the kernel algebra were constructed 
for the class ([T]) in [l8j in the course of partial preliminary group classification involving a ten- 
dimensional subalgebra of the equivalence algebra of this class. All these cases are G~-equivalent 
to particular subcases of Cases 1-5 from Table [1] of the present paper. 

The approach of complete preliminary group classification can be optimized via selecting of 
appropriate subalgebras of i.e. subalgebras whose projections to the space of system vari- 
ables are maximal Lie invariance algebras of systems from the class C\s- The simplest common 
property of appropriate subalgebras is that they contain the kernel algebra. Other criteria for 
selecting of appropriate subalgebras including bounds for dimensions of extensions or additional 
extensions are derived via examination of the determining equations for Lie symmetries of sys- 
tems from the class C\s- In a certain sense, this means combining the algebraic method of group 
classification with the direct method based on the study of compatibility and the integration 
of the determining equations up to G~-equivalence. The usage of the optimized technique of- 
ten allows one to reduce the classification problem to classification of certain low-dimensional 
subalgebras of the equivalence algebra, even if the equivalence algebra is infinite-dimensional 
and there exist infinite-dimensional extensions of the kernel. Related calculations are not too 
cumbersome. Thus, minimal computations which are necessary for complete preliminary group 
classification of class ([1]) are exhausted by the first parts of Sections H] and [8] and entire Sec- 
tions [Sj El [10] and [TTJ These computations result in the absolute majority of inequivalent cases 
of Lie symmetry extensions for the class ([1]) , which are presented in Table [1] (the exceptions are 
only Cases 6, 18, 19 and 22). 

There exist two ways to apply the algebraic method to complete group classification. The 
first way is to reduce complete group classification to preliminary one. The reduction can be 
realized, e.g., via proving that the class C\s is weakly normalized or partitioning this class into 
weakly normalized subclasses and other subclasses which can be easily classified by the direct 
method. Although the partition into subclasses usually involves cumbersome and sophisticated 
computations, it is an effective tool of group analysis since it accurately adapts the classification 
procedure to the structure of the class C\s- This is the way that has been used in the present 
paper. The class ([I]) is partitioned into two subclasses possessing the same equivalence group as 
the whole class ([1]). One of the subclasses is normalized, the other is semi-normalized and mapped 
by equivalence transformations onto its subclass (|29]) of simple structure. Group classification 
of the subclass (|29p has been obtained in the course of the partition which results in only four 
special cases of Lie symmetry extension (Cases 6, 18, 19 and 22 from Tabled]), which are not 
related to subalgebras of £|~. The second way is to directly classify G~-inequivalent appropriate 
algebras contained in the span g^^ = (g^l^ € S) of maximal Lie invariance algebras, Qq, of 
all systems from the class C.\s- This way properly works only if the class C\s possesses certain 
properties, e.g., if the maximal Lie invariance algebra Qg is of low dimension for any 6 ^ S [22J or 
if the class C\s is at least weakly normalized or partitioned into weakly normalized subclasses but 
this property is not explicitly checked [3] 127 1 I^Sl [54j . An explanation for the above observation 
is that G~-equivalence is not appropriate in the course of classification of subalgebras contained 
in g^-* if g^^ is strongly inconsistent with the equivalence algebra g~ (e.g., much wider than the 
projection Pg'" of g~ to the space of system variables). 

Due to the above partition of the class ([1]) we have obtained essentially stronger results than 
the solution of the usual group classification problem by Lie-Ovsiannikov for this class. The 
partition exhaustively describes the set of admissible transformations in the class ([1]). Moreover, 
the fact that the whole class ([I]) is semi-normalized guarantees that there are no additional point 
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equivalence transformations between cases of Lie symmetry extensions presented in Table [H i.e., 
the same table gives the complete group classification of the class ([1]) with respect to general 
point equivalence. 

The extension and clarification of the group classification toolbox is by no means a pure 
mathematical problem. Methods from symmetry analysis including group classification have 
the potential to provide, in particular, a unifying framework to construct invariant local closure 
or parameterization schemes for averaged nonlinear differential equations \33\ 140} |47] . As finding 
appropriate closure ansatzes for averaged differential equations is at the basis of any numerical 
model of (geophysical) fluid dynamical systems, it is immediately clear that group classification 
can play a crucial role in the construction of different computational codes for such systems. 
The classes of differential equations arising in the course of the parameterization problem are 
usually much wider and have more complicated structure than the classes studied in conventional 
group classification. It generally cannot be expected to completely solve the group classification 
problems for such classes using existing methods. Hence the development of new tools for 
group classification of differential equations simultaneously with the improvement of well-known 
approaches still remains an attractive and challenging problem. Especially for the above complex 
classification problems, the whole framework of the algebraic method as described and extended 
in the present paper including the proposed algebraic method of finding equivalence groups 
seems to be most appealing. 
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